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20.  ABSTRACT  (ContlniMd) 

prtp«r«t1on  for  futurt  lnifrorwwnts  In  WBMOO,  Analysis  and  intarprata- 
tlon  of  the  Wideband  Satellite  data  base  from  Poker  Flat,  Alaska  was  continued. 
The  utility  of  the  phase-gradient  ratio  (l.a.,  the  ratio  of  nas  phase  difference 
for  InterferoMter  baselines  perpendicular  and  parallel  to  known  sheetllke 
Irregularities)  was  Investigated  for  detenalnlng  Irregularity  axial  ratios  and 
was  found  to  be  of  potential  use.  The  ratio  varies  strongly,  however,  with  the 
phase  spectral  liwlex,  p.  Thus,  quantitative  use  of  the  phase-gradient  ratio 
depends  upon  prior  evaluation  of  the  geophysical  and  geoantrlcal  behavior  of  p.C 
A  study  was  Initiated  of  the  b^avlor  of  p  In  the  data  base  from  Poker 
Flat.  The  results  show  a  Marked  variation  with  the  angles  between  the  line  of 
sight  and  the  local  L  shell  and  Magnetic  Meridian.  This  behavior  Is  Intei'preted 
tentatively  as  evidence  for  scale>s1ze  dependence  of  the  degree  of  cross-field 
anisotropy  exhibited  by  the  Irregularities,  tending  froM  sheetllke  structures  at 
the  largest  scales  to  rodllke  structures  at  the  SMallest.  In  the  course  of  the 
phase  spectral  study,  two  persistent  features  were  found.  The  first  Is  a  downward 
break  In  thw  high-frequency  spectral  tall  at  about  10  Hz  (»  300  meters  scale  size 
In  the  direction  of  scan,  which  Is  roughly  north-south  geomagnetic).  This 
feature,  which  was  noted  In  about  a  quarter  of  the  spectra  Inspected,  Is  similar 
to  breaks  reported  In  equatorial  data.  The  second  feature,  which  occurs  about 
three-quarters  of  the  time.  Is  an  upward  break  on  the  low-frequency  end  of  the 
spectrum  (below  0.5  Hz,  which  corresponds  to  about  6  km). 
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SECTION  1 


INTRODUCTION 

In  carrying  out  Its  mission  to  evaluate  a  wide  variety  cf  nuclear  effects  on  weapons 
systems,  the  Defense  Nuclear  Agency  (ONA)  conducts  a  vigorous  research  program  on  h  igh-altltude 
plasmas  and  their  Impact  on  radio  propagation.  In  recent  years,  this  program  has  concentrated 
on  medium-scale  structures  (Keskinen  and  Ossakow,  1902)  and  the  potentially  disruptive  radio¬ 
wave  scintillations  that  they  produce  (Wittwer,  1979).  A  major  element  of  that  program's  focus 
was  DNA's  Wideband  Satellite  Experiment  (R^no  et  ol,  1977;  Fremouw  et  al,  1978). 

High-latitude  results  from  Wideband  were  put  Into  a  form  directly  useable  for  engineering 
evaluation  of  auroral  zone  scintillation  effects  by  means  of  a  plasma-irregularity  model  and 
propagation  routines  encoded  In  a  computer  program  called  WBMOO  (Fremouw  and  Lansinger,  1981). 
The  Wideband  data  base  from  Poker  Flat,  Alaska,  was  used  for  Uerative  development  and  testing 
of  the  model-  The  prograiii,  written  In  Fortran,  was  constructed  for  highly  Interactive  applica¬ 
tion  by  a  systems -oriented  user. 

This  report  summarizes  recent  Improvements  In  WBMOO.  It  also  describes  research  Into 
Irregularity  characteristics  revealed  by  the  Poker  Flat  data  base,  which  could  lead  to  further 
Improvement  In  the  high-latitude  model  and  to  guidance  for  DNA's  forthcoming  KILAT  satellite 
program  (Fremouw,  1982). 

The  code  Improvements  are  summarized  In  Section  2,  starting  with  a  routine  to  order  hlgh- 
latltude  Irregularity  characteristics  in  geomagnetic  rather  than  geographic  time.  Correction 
of  a  subtle  error  in  treating  the  llne-of-sight  scan  velocity  of  an  orbiting  satellite  In  the 
phase-screen  propagation  theory  employed  Is  described  In  Section  2-2,  followed  by  refinement  of 
the  orbit  calculation.  Section  2-4  summarizes  Implementation  of  WBMOO  for  operational  use  at 
USAF  Global  Weather  Central,  as  fully  documented  elsewhere  (Secan,  1982).  A  major  Improvement 
for  application  of  WBMOO  to  use  of  geostationary  satellite  links  at  high-latitude  stations  Is 
described  In  Section  2-5,  namely.  Incorporation  of  a  realistic  model  (Heells  et  at,  1982)  for 
convective  drift  of  Irregularities  In  the  high-latitude  F  layer. 

A  three-character  revision  numbering  system  has  been  adopted  to  Identify  major  changes  to 
the  WBMOO  program.  The  first  (numeric)  character  Identifies  revisions  In  the  overall  structure 
of  the  code  that  Is  not  directly  part  of  the  Ionospheric  Irregularity  model.  The  second 
(alphabetic)  character  Identifies  changes  In  the  algebraic  form  of  the  Irregularity  model 
beyond  simple  revisions  In  the  numerical  constants  that  calibrate  the  model.  The  final 
(numeric)  character  Identifies  changes  in  these  numerical  constants.  Lesser  changes  will  be 
noted  only  by  changes  In  the  revision  date. 

This  three-character  revision  number  appears  In  a  comment  near  the  beginning  of  Program 
WBMOO,  In  the  Interactive  dla.ogue  with  the  user,  and  in  the  main  output  file.  The  ntodel 
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contained  in  the  code  (e.g.  B3)  will  be  identified  In  comeent  statements  near  the  beginning  of 
rvibroutine  MDLPRM  and  function  FCSL.  The  current  version  of  UBHOD,  Revision  6B3,  Is  summarized 
In  Section  3. 

Research  topics  are  discussed  In  Section  4,  starting  with  exploration  of  efficient  means 
for  deducing  the  three-dimensional  configuration  of  scinti 11  at Ion -producing  irregularities 
from  Interferometer  data.  One  such  t'zchnlque  was  found  to  require  Independent  Information  on 
the  one-dimensional  spatial  spectrum  of  the  Irregularities,  and  spectral  behavior  Is  discussed 
In  Sections  4-2  and  4-3. 
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SECTION  2 


INPROVEMENTS  IN  WBHOD 


2-1  MAGNETIC  TIME. 

Most  high-latitude  Ionospheric  phenomena  are  better  ordered  diurnal 1y  In  magnetic  time 
than  In  geographic-meridian  time.  Accordingly*  the  expression  coded  In  WBMOO  to  describe  the 
magnetic  Invariant  latitude,  Xjj,  of  the  high-latitude  scintillation  boundary  Is  as  follows: 


Xu  •  71.8“  -  1.5K, 


Tr(t  -  2) 
-  5.5  cos  — — 


(1) 


where  Kp  Is  the  planetary  magnetic  activity  Index,  and  t^  Is  time  In  hours  past  magnetic 
midnight. 

For  Iterative  testing  against  data  from  Poker  Flat,  a  fixed  offset  was  used  between 
magnetic  and  geographic  meridian  time  (namely,  magnetic  midnight  lagging  geographic  meridian 
midnight  by  1.5  hours).  For  use  of  the  code  at  other  locations.  It  Is  necessary  to  establish 
the  magnet 1c -midnight  offset  as  a  function  of  geographic  coordinates.  A  simple  routine  for 
calculating  this  offset  was  developed  and  Implemented  in  Rev  3A2  and  later  modified  In  Rev  6B3. 
Both  methods  employed  will  be  presented  here. 

2-1.1  Revision  3A2. 

The  basis  for  the  magnetic-time  calculation  Is  the  set  of  definitions  established  by 
McNIsh  (1936).  Magnetic  midnight  Is  defined  as  the  time  at  which  a  line  from  the  sun  through 
the  earth's  center  Intersects  the  night  side  of  the  earth  on  the  station's  geomagnetic  merid¬ 
ian.  Magnetic  time  Is  then  taken  to  advance  at  a  constant  rate  of  24  hours  per  revolution  of 
the  earth  on  Its  spin  axis.  For  calculating  the  magnetic  longitude,  of  a  station,  the 
following  dipole-based  expression  Is  used: 

tan  ”  -tan(t  -  i)  sin  X  sec  (x  +  X  )  (2) 

In  0  0 

where  x  “  tan"^[cos()l  -  !Q)ctn  x] 

e.,X  >  geographic  longitude  and  latitude  of  station 

and  tp,  Xq,  -  geographic  longitude  and  latitude  of  geomagnetic  poles. 

The  calculation  of  the  magnetic-time  offset  begins  with  a  call  by  Subroutine  CCGLT  to 
Subroutine  CGLT  for  calculation  of  the  geomagnetic  longitude  of  an  ionospheric  penetration 
point,  using  Eq.  (2).  A  second  call  to  CGLT  finds  the  geomagnetic  longitude  of  -he  effective 
anti-sun  at  the  time  that  it  lies  In  the  penetration  point's  geographic-meridian  plane.  The 
anti-sun  geomagnetic  longitude  Is  then  adjusted,  by  means  of  an  Iterative  loop,  to  match  that 
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of  the  penetration  point.  The  corresponding  adjusbnent  In  geographic  longitude,  expressed  1n 
hours  at  the  rate  of  15®  per  hour,  constitutes  the  lead  time  of  magnetic  midnight  relative  to 
geographic  meridian  midnight. 

To  calculate  the  true  geomagnetic  longitude  of  the  physical  sun,  one  would  enter  Subrou¬ 
tine  CGLT  using  the  decllnetlon  of  the  sun  for  the  Input  geographic  latitude.  There  would  then 
be  a  further  correction  for  the  equation  of  time.  To  calculate  the  true  geomagnetic  longitude 
of  the  anti-sun,  one  would  employ  the  negative  of  solar  declination  for  the  Input  latitude  and 
then  make  the  Identical  equatlon-of-time  correction. 

Using  the  solar  declination  and  correcting  for  the  equation  of  time  both  Introduce 
seasonal  variations  In  the  geomagnetic  longitude  of  the  sun  for  a  given  geographic  longitude. 
The  variation  Introduced  by  the  equation  of  time  Is  the  same  for  the  sun  and  the  anti-sun.  It 
Is  easy  to  show,  however,  that  the  seasonal  variation  Introduced  by  employing  the  true  declina¬ 
tion  Is  opposite  for  the  sun  and  anti-sun.  It  Is  curious  but  true,  therefore,  that  the  seasonal 
variation  In  magnetic  time  as  commonly  defined  depends  upon  whether  one  chooses  noon  (the  sun) 
or  midnight  (the  anti-sun)  as  the  time  origin. 

One  may  avoid  the  difference  Imposed  by  the  above  arbitrary  choice  by  employing  zero 
declination  for  either  the  sun  or  the  anti-sun,  which  corresponds  to  equinoctial  conditions, 
and  we  have  done  so.  To  do  otherwise  would  be  to  Introduce  an  Ill-based  seasonal  variation  Into 
our  scintillation  model.  Since  we  hope  to  Identify  and  code  a  true  seasonal /longitudinal 
variation  1n  scintillation.  Introduction  of  such  a  variation  a  priori  (and,  In  effect,  inadver¬ 
tently)  would  Introduce  confusion  Into  our  efforts.  For  this  reason,  we  have  coded  zero  as  the 
declination  of  the  anti-sun,  and  we  have  chosen  also  to  Ignore  the  seasonal  correction  due  to 
the  equation  of  time.  We  shall  keep  the  behaviors  described  thereby  In  mind  as  we  investigate 
seasonal/longitudinal  variation. 

Figure  1  illustrates  the  difference  between  magnetic  and  geographic-meridian  times  com¬ 
puted  by  means  of  the  magnetic-time  offset  routine.  Magnetic-time  advance  Is  presented  as  a 
function  of  geographic  longitude  for  four  geographic  latitudes. 

2-1.2  Revision  6B3. 

The  procedure  to  convert  local  mean  solar  time  (t^)  to  local  magnetic  (dipole)  time  (t^) 
was  modified  In  Rev  6B3  to  allow  direct  rather  than  Iterative  calculation  of  t^.  First,  the 
dipole  longitude  Is  calculated  from 

cosXs1n(J,- 1^) 

^m  ”  sinx^cos  cos(t- t^j-cos  X^sln  X 

where  all  variables  are  as  defined  In  Equation  2.  The  geographic  longitude  { at  which 
crosses  the  geographic  equator  Is  calculated  from 

“  ^0  tan  ij.  (4) 
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Geographic 

Latitude 


GEOGRAPHIC  LONGITUDE  •  deg  east 


Figure  1.  Difference  between  geomagnetic  time  and  ge*,graph 1c- meridian 
time  as  function  of  geographic  latitude  and  longitude.  Neg¬ 
ative  values  mean  that  geomagnetic  midnight  lags  geographic 
meridian  midnight. 


The  correction  to  LT  then  Is  defined  as  the  difference  between  l  and  i  , 
£  -  t  e 

t  »  t  +  — 
m  '’s 


This  new  algorithm.  Implemented  In  Function  GMLTC,  yields  the  same  results  as  subroutines  CCGLT 
and  CGLT  described  In  Section  2-1.1,  but  It  Is  quicker,  as  no  Iteration  Is  required. 

An  additional  refinement  to  this  calculation  is  use  of  the  dipole  pole  location  (X  ,9.) 
calculated  from  lGRF-1980  using  the  expressions 


sin 
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■  9l0  *  ’ll  *  •'ll 

tan  i  =  ^11 
9ll 

where  and  hp,^  are  the  (Schmidt  quasi -normalized)  spherical  harmonic  coefficients  for  the 
ISRF  expansion.  For  epoch  1980.0  these  values  are  =«  78.8®N,  =  70.76"W. 

2-2  SCAN  VELOCITY 

During  the  course  of  iteratively  testing  WBMOD  Rev  4A2  against  Wideband  scintillation  data 
from  Goose  Bay,  we  encountered  a  problem  in  calculating  scan  velocity  of  the  line  of  '■''ht 
between  a  low-orbiting  satellite  and  a  fixed  ground  station.  The  problem  seriously  affected 
scintillation  calculations  for  the  observing  geometry  at  Goose  Bay,  but  had  not  been  noticed 
for  the  Poker  Flat  geometry.  The  problem  was  traced  not  to  a  program  bug  but  to  a  subtle  error 
in  the  treatment  of  scan  velocity  in  the  scattering  theory  employed  (Rino,  1979). 

Prior  to  discovery  of  the  problem,  velocity  was  calculated  in  a  subroutine  called  VXYZ 
using  a  coordinate  system  centered  at  the  ionospheric  penetration  point  (Rino  and  Fremouw, 
1977).  The  geometry  was  developed  from  the  point  of  view  of  spaced-receiver  measurements  on  an 
observing  plane.  For  application  to  single-receiver  measurements  of  scintillations  in  the  time 
domain,  the  approach  was  to  calculate  the  velocity  of  the  receiver  as  seen  from  the  location  of 
the  ionospheric  penetration  point.  The  problem  was  that  although  the  theory  accounted  for 
translation  of  the  coordinate  origin  relative  to  the  receiver,  it  did  not  take  into  account  a 
rotation  of  the  coordinate  system  relative  to  the  local  vertical. 

For  an  observer  located  at  the  ground-based  receiver,  the  moving  coordinate  system, 
coupled  with  the  fact  that  irregularity  anisotropy  is  handled  in  the  scattering  theory  by  means 
of  an  effective  velocity,  Vg,  defined  in  a  nonisotropic  coordinate  system,  makes  the  velocity 
calculation  somewhat  obscure.  The  situation  is  quite  clear,  however,  from  the  observer's  point 
of  view  in  a  fixed  reference  frame  centered  at  the  receiver,  especially  for  isotropic  irregu¬ 
larities.  The  relevant  velocity  then  obviously  is  that  of  the  line  of  sight  perpendicular  to 
itself  at  the  distance  of  the  penetration  point.  Thus,  for  isotropic  irregularities,  Vg  must 
reduce  to 

Vj^  =  cos^e  (6) 

where  R  is  the  range  from  the  receiver  to  the  penetration  point,  e  is  the  elevation  angle  of  the 
source  as  seen  from  the  receiver,  and  a  is  the  source  azimuth  from  the  receiver.  We  found  that, 
for  the  Goose  Bay  observing  geometry,  it  did  not,  as  illustrated  in  Figure  2. 
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Figure  2.  Comparison  of  Ve  (dashed)  calculated  In  WBMOD  Rev  4A2  for  Isotropic 
Irregularities,  with  Vj^  (solid)  calculated  In  receiver-centered 
coordinate  system  (Equation  1).  The  former  should  equal  the  latter, 
but  does  not.  It  Is.  therefore.  Incorrect. 


The  difficulty  stems  from  rotation  of  the  coordinate  system;  It  may  be  appreciated  from 
Figure  3,  which  Illustrates  the  geometry  at  a  calculation  time,  t^,  and  at  two  time  Increments, 
±At.  The  scattering  theory  Invokes  a  locally  plane  phase  screen  oriented  horizontally  at  the 
tg  penetration  point.  The  geometry  Is  calculated,  however,  from  a  coordinate  origin  moving 
along  a  surface  of  constant  height  above  a  spherical  earth. 

The  theory  requires  a  two-dimensinal  velocity  referred  to  the  Instantaneous  (planar)  phase 
screen,  which  Is  given  by 

Vjjj  "  ® 

Vjy  •  Vy  -  tan  9  sin  «(•  (7b) 

where  9  and  are  the  Incidence  angle  and  magnetic  heading,  respectively,  of  the  propagation 

vector,  and  V  ,  V  ,  V  are  the  three-dimensional  velocity  components  of  the  receiver  In  the 
A  y  z 

po1nt«centered  coordinate  system.  Equations  (7)  correctly  shift  the  origin  from  the  curved  z 
■  0  surface  to  the  Intersection  of  the  rays  with  the  planar  phase  screen.  Ihey  do  not  yield  the 
correct  two-dimensional  velocity,  however,  unless  the  coordinate  rotation  Is  accounted  for 
Independently. 
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Th«  coordinate  rotation  Is  Indicated  In  Figure  3  by  the  small  angles,  end  The 
velocity  error  arose  from  calculating  the  receiver  position  in  a  system  oriented  along  the 
local  vertical  at  the  ±  points.  Indicated  by  primed  coordinate  distances  and  the  dashed 
triangles  In  Figure  3.  The  correct  velocity  Is  obtained  by  calculating  the  receiver  position 
In  a  system  oriented  parallel  to  the  Instantaneous  vertical  at  the  calculation  time,  t^,  as 
Indicated  by  the  nonprimed  coordinate  distances  and  the  solid  triangles.  Our  finite- 
differences  solution  Mas  to  rotate  through  the  angle,  l  ,  from  the  primed  to  the  non-primed 
coordinate  system,  by  means  of  the  fol lowing  equations: 


■  (x^  sin  ♦  2^  cos  Xj)  cos(i!.(,-  -  yj  s1n(tjj-  sin  X^ 

+  (x^  cos  x^  -  2^  sin  X^)  cos  X^  (8a) 

y^  -  (x^  sin  X^  +  2^  cos  >^)  slnCt^-t  )  ♦  y^  cos(t^-£j  (8b) 

2^  ■  (x^  sin  X^  ♦  2^  cos  X^)  cos(«.^-t^)  -  y^  sln(t^-t^)  cos  X^ 

-  (Xj  cos  X,.  -  2^  sin  X^)  sin  X^  (8c) 


where  the  primed  and  unprimed  coordinates  are  as  In  Figure  3,  and  x  and  t  respectively 
represent  (geographic)  latitude  and  longitude  of  the  penetration  points  Indicated  by  the 
subscripts. 

In  the  case  of  spaced -receiver  measurements.  Equations  (7)  would  account  for  the  tilt 
angle  of  the  observing  plane  (remote  section  of  the  earth's  surfr.ce)  In  the  point -centered 
cocrilnate  system.  In  the  situation  at  hand,  it  accounts  only  for  the  small  angles  between  the 
2-0  surface  and  the  Instantaneous  planar  phase  screen,  which  in  fact  turns  out  to  be 
negligible.  The  other  effect  of  surface  curvature  (coordinate  rotation)  produced  substantial 
velocity  errors  In  certain  circumstances,  and  It  Is  now  corrected  for  by  means  of  Equations 
(8). 

The  In-screen  velocity,  Vj,  Is  required  because  the  scattering  theory  effectively  col¬ 
lapses  the  thick  Irregularity  layer  as  a  shadow  pattern  onto  the  equivalent  thin-phase  screen. 
The  relevant  physical  velocity,  however,  still  Is  that  perpendicular  to  the  line  of  sight.  The 
angular  component  of  velocity  Is  extracted  and  the  range-rate  component  suppressed  by  means  of 
the  following  expression  for  the  effective  velocity,  V^: 


(CV 


sx 


•  »Vs,  * 

(AC  -  b2/4)1/2 


(9) 
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Anisotropy  Is  accounted  for  (Rino  and  KreMouw,  197>)  by  means  of  A,  3,  and  C,  which  are 
functions  of  0  and  ^  and  of  the  Irregularity  axial  ratios,  a  and  b,  and  orientation  angles 
(magnetic  dip  and  off-L-shell  angle). 

Figure  4  shows  calculated  from  Equation  (9)  In  WBHOO,  employing  the  new  velocity 
subroutine.  The  solid  curve  Includes  the  effects  of  an1sotrt)py  Inherent  In  the  Irregularity 
model,  and  the  dashed  one  was  calculated  with  a  and  b  set  to  unity.  Comparison  of  the  latter 
with  the  solid  curve  In  Figure  2  shows  that  It  Is  Identical  to  calculated  from  Equation  (6), 
as  It  should  be.  We  also  calculated  the  nonIsotropIc  V^,  In  the  receiver-centered  coordlna^ie 
system,  employing  only  the  angular  components  of  scan  velocity,  which  Is  Intuitively  the  clear 
approach.  The  result  Is  Identical  to  the  solid  curve  In  Figure  4,  indicating  that  the  range- 
rate  velocity  component  has  been  properly  suppressed  In  the  latter  in  the  presence  of  aniso¬ 
tropy. 

In  addition  to  results  from  the  new  velocity  subroutine.  Figure  4  shows  (dotted  curve) 
calculated  In  WBHOO  prior  to  Implementing  the  velocity  correction  described  In  the  foregoing. 
The  resulting  curve  approximates  the  correct  one  (solid)  through  much  of  the  Goose  Bay  satel¬ 
lite  pass  employed,  but  It  departs  substantially  at  the  higher  penetration-point  Invariant 
latitudes.  The  high-latitude  dip  In  effective  velocity,  which  caused  us  to  suspect  a  problem 
in  the  first  place,  was  Indeed  the  result  of  an  error. 


I 

8 

_J 


Figure  4.  calculated  from  new  velocity  subroutine  for  anisotropic  (solid)  and 
Isotropic  (dashed)  Irregularities.  Note  Identity  of  latter  to  solid 
curve  In  Figure  2.  Dotted  curve  represents  Incorrect  Vg  calculated  for 
nonIsotropIc  Irregularities  uslno  old  velocity  subroutine. 
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The  cnnfMiutnc*  of  tho  velocity  error  on  coeputatlon  of  phase  scintillation  was  sub¬ 
stantial  In  the  Goose  B«y  pass  at  han<ll.  The  difference  between  the  phase  scintillation  Indices 
calculated  fron  the  saw  Irregularity  aodel  using  the  nsw  and  old  velocity  routines  Is 
strikingly  shown  In  Figure  5.  The  totally  unexpected  and  peculiar  null  In  nas  phase  fluctua¬ 
tion  calculated  with  the  old  routine  resulted  solely  fron  the  velocity  error. 

The  effect  of  the  velocity  error  Is  quite  different  for  different  pass  geowtrles, 
probably  having  close  to  Its  greatest  lepact  In  the  geowtry  underlying  Figures  2,  4,  and  5.  It 
Is  Much  less,  for  Instance,  In  Figure  6,  which  contains  phase  scintillation  Indices  calculated 
using  the  old  (dotted)  and  new  (solid)  velocity  subroutines  for  a  nIghttiMe  Wideband  pass 
essentially  along  the  Magnetic  wridlan  over  Goose  Bay.  The  discontinuity  In  the  dotted  curve 
stewed  fron  Insufficient  nuwrical  accuracy  In  the  old  velocity  subroutine  as  the  penetration 
point  Moved  across  the  Magnetic  wridlan.  It  too  was  renedled. 


Figure  5.  Phase  scintillation  Index  (ms  fluctuation)  calculated  In  WBHOD,  using  old 
(dotted)  and  new  (solid)  velocity  subroutines,  for  nighttime  pass  of  Wide¬ 
band  satellite  east  of  Goose  Bay. 
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INVARIANT  LATITUDE  (deg) 


Figur.  6.  Ftg^  "rti"'  •*“  “•  ■•S'**’'  ■•'<'»•'’ 

■j  «®ose  Bey.  Agetn«  tee  solid  end  dotted  curves  resoectlvulv  sKom  raioiia. 
tions  using  the  nev  end  old  velocity  subroutines.  celcOe- 


2*3  ORBIT  CALCULATION 

During  entlysis  of  the  seen  velocity  proble*  discussed  In  Section  2-2  end  concurrent  work 
on  definition  of  the  orbit  of  the  HILAT  setelllte,  elgorlthn  were  developed  thet  would  el low 
direct  celculetlon  of  e  circuler  orbit  frow  two  Input  setelllte  locutions.  These  elgorlthes, 
lepleMnted  In  subroutines  STORB  end  LATLON.  heve  repleced  subroutines  FN00R8,  SRCHH,  end  SRCKT 
end  function?  SUT  end  SLON  (Rino  et  al,  1978). 

Subroutine  STORB  celculetes  the  necessery  orbitel  pereeeters  fro*  input  velues  of  setel¬ 
llte  eltitude  (hj).  two  setelllte  locutions  (X^.  t  and  X-.  ^),  end  the  time  et  the  first 
locution  (tj)  es  follows: 

1.  Celculete  the  orbitel  enguler  velocity  of  the  setelllte  from 
n  •  3600  -It  redlens/hour 

S 

where 

Rj  •  hs  +  6371.2  km 
u  -  3.986013x10®  km®/sec^. 

2.  Iteretively  celculete  the  time  et  (x^.ij)  from  the  orbitel  enguler  velocity,  the 

enguler  distence  between  (Xj.lj)  end  end  the  eerth's  rotetlonel  enguler  velocity  (0  • 

^redlens/hour).  * 


3.  Correct  the  second  longltode  for  the  earth's  rotation  2  ”  z  *  e^^2  ~ 
calculate  the  longitude  (  and  tint  (t^)  of  the  last  ascending  node  and  the  orbital  Incllna* 
tion  angle  (1)  In  a  non-rotating  frane  as  follows: 


l 

0 


I 

tan  Xj,  s1n(t2-<^X^ 
Un  Xj  -  tan  tjCos(t2-tj) 


(10) 


s1n( 

£in  ^ 


t  .  t  -  ^ 

^0  ^1  n 

where  4>  Is  the  angular  d1 


I 

stance  along  the  orbit  fron  (0,  t^)  to  (X^t^). 


(11) 

(12) 


4.  Revert  to  the  rotating  frane  by  correcting  with 

t  .  t'  -  • 

0  0  =  «  * 

The  orbital  paraneters,  n,  1,  and  t^,  are  used  by  subroutine  LATLON  to  calculate  the  sub- 
satellite  location  (X,<!)  for  a  given  t1ne»  t,  as  follows: 


1.  Calculate  the  angular  distance  (♦)  along  the  orbit  from  (O.t^)  using 
$  ■  n(t-t^)  . 

2.  Calculate  the  location  fron 

X ■  s1ir^(s1n  1  sin  ♦ )  (13) 


Hq  +  tan 


-1 


cos  1  sin  ♦ 
cos 


(14) 


These  algorithms  allow  the  user  to  generate  a  locally  circular  orbit  that  will  approximate 
the  trajectory  of  any  actual  satellite  over  small  suction  of  Its  orbit. 

2-4  IMPLEMENTATION  AT  AFGWC 

One  of  the  tasks  carried  out  was  to  modify  the  most  recent  version  of  WBMOO  (Rev  6A3)  for 
operational  use  by  the  USAF  Air  Heather  Service  (AHS)  at  the  Air  Force  Global  Heather  Central 
(AFGHC).  Completion  of  this  task  required  the  following  work: 

1.  Meet  with  AHS  representatives  at  AFGHC  to  determine  the  requirements  to  be  met  by  an 
operational  version  of  HBMOO.  These  requirements  were  formalized  In  a  Functional  Description 
(Secan,  1982)  and  a  project  development  plan. 

2.  Convert  the  existing  HBMOD  code  to  meet  the  ANSI  X3. 9-1978  FORTRAN  standard  and 
Insert  standardized  In-line  documentation. 
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3.  Otvelop  AF6MC-sptc1flc  1nt«rfac«  routines  to  Interact  with  the  AFGHC  Spacn  Environ* 
litnt  Support  Branch  (AFSUC/USE)  Operations  Center  CRT  and  Data  CoMMunIcatlons  Terminal  (DCT) 
and  the  AFGWC  Astrogeophysical  Data  Base  (AGDB).  Figure  7  shows  an  example  of  the  Input 
Interaction  between  the  AF6UC  version  of  WBNOO  and  the  user.  Figure  8  Is  a  sample  of  the 
output,  which  can  be  routed  by  the  user  tc  a  local  off*11ne  printer  (DCT). 

4.  Establish  procedures  for  Physical  Dynamics  to  provide  AFGWC  with  model  updates  as 
they  occur. 

5.  Write  and  publish  system  documentation  In  accordance  with  DoD  Standard  7935. 1-S 
(1977)  (Secan.  1982). 

6.  Implement  and  test  the  reconfigured  WBHOD  at  AFGWC.  This  Included  on-site  training 
of  AFGWC  operations  and  programming  staff  members. 

This  task  was  begun  on  November  17,  1981,  with  a  meeting  at  AFGWC  and  was  completed  on 
August  17,  1932,  when  the  following  Items  were  provided  to  AFGWC: 

1.  Documentation  required  by  DoD  Standard  7935. 1-S  (Secan,  1982) 

Vol.  I.  Functional  Descriptio-; 

Vol.  II.  Users  Nanusl 

Vol.  Ill  Program  Maintenance  Manual 

Vol.  IV  Test  Plan 

Vol.  V  Test  Analysis  Report 

?.  A  magnetic  tape  containing  the  latest  version  of  WBMOD  (6A3)  as  configured  for  use  at 
AFGWC. 

3.  A  list  of  the  tape  contents. 

4.  Basell.ie  output  samples  for  the  non-equatov ial  cases  listed  In  Appendix  A  of  the 
Program  M,i1:itenance  Manual. 

All  sections  o^  W8M00-6A3  were  Iwpimnented  at  AFGWC,  with  the  exception  of  equatorial -zone 
scintillation.  This  section  of  the  model  required  further  calibration.  Rather  than  to  provide 
scintillation  estimates  that  were  unreliable,  1i  was  decided  to  "switch  off"  the  model  inter¬ 
nally  within  25  '  of  \,he  dip  equ't&r.  All  Involved  personnel  at  AFGWC  are  aware  r*  this 
temporary  constraint. 

2-5  IRREGULARITY  DRIFT  VELOCITY 

2-5.1  Introduction 

The  WBMOD  program  has  developed  to  where  It  provides  reliable  estimates  of  average  scin¬ 
tillation  strength  tor  trans Ionospheric  radiowave  systems  with  one  terminus  In  low  earth  orbit. 
It  contained  a  serious  deficiency,  however,  for  application  to  systems  having  a  very  high 
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IONOSPHERIC  SCINTILLATION  HODEL 
(VERSION  6A3  -  18  FEB  1982  ) 

ENTER  initial  UT  DATE  (NH  DO  YY> 

11  29  80 

ENTER  FIO  VALin  FOR  29  NOV  80 
190 

ENTER  INITIAL  UT  TIME  (HHMHl 
0000 

««AP  IS  NOT  AVAILABLE  IN  THE  AODB*« 

ENTER  AP  VALID  FOR  29/03O0UT  NOV  60 

40 

ENTER  SYSTEM  OPERATING  FREOUENCY  (MHZ) 

137.68 

ENTER  PHASE  STABILITY  DURATION  TIME  (SEC) 

10 

ENTER  receiver  LOCATION 

LAT(-H«RTH)  L0N(4EAST)  ALTITUOE(KM) 

69.  1  >147.  9  0.  2 

ENTER  INITIAL  SATELLITE  POSITION 
LAT(4NQRTH)  L0N(’*EAST>  ALTITUDE(KM) 

0.0  -80.0  42240.0 
ENTER  RUN  MODE  DESIRED 
1:  ORBIT  2:  STEP  PARAMETER  3:  STEP  RECEIVER 

2 

ENTER  PARAMETER  TO  VARY 

l:TIME  2:AP  3; FIO  4:  FREQ  9;  SATELLITE  LONGITUDE 
9 

ENTER  FINAL  SATELLITE  LONGITUDE  (•►EAST} 

-180.  0 

ENTER  NUMBER  OP  INCREMENTS  (90  MAX) 

11 

DO  YOU  WANT  A  SUMMARY  OF  INPUTS?  (1:YE8.  2:  NO) 

1 

OPTION:  STEP  PARAMETER 
DATE:  11/29/80  TIME:  OOOO  UT 
FIO:  190  AP:  40 
SYSTEM  PARAMETERS 

FREOUENCY:  137.7  MHZ 

RECEIVER  LOCATION:  69.  t/-147.  9/  0.  20  KM 
satellite  LOCATION:  0.  0/  -80.  0/42240.  O  KM 

ONE-WAY  PROPAGATION 
PHASE  STABILITY  PERIOD:  10  SEC 
VARIED  PARAMETER  (SLON) 

FINAL  VALUE:  -180.  GO 
NUMBER  OF  STEPS:  11 

DO  YOU  WANT  TO  CHANGE  ANY  1)«>UT3?  (1:  YES.  2:  NO) 

2 

[Output  List  -  S««  Figure  8] 

DO  YOU  WANT  A  DCT  COPY?  (1:YES.  2:  NO) 

2 

DO  YOU  WANT  ANOTHER  RUN?  (1:YE8.  2:  NO) 

2 

WBMOD: FIN 


Figure  7.  Sample  Input  Interaction  (Orbit  Mode) 
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lONc 

3PHERIC  SCINTILLATION  MODl£L 

(VERSION 

6A3  -  18 

FEB  1982  ) 

DAiIi 

:  06/13/62 

TIME: 

2230  FIO: 

130  AP: 

40  FREQUENCY: 

137.  7  MHZ 

RECEIVER; 

69.  IN.  147.  9U.  0.  SOKM 

SATELLITE:  80.  ON.  89.  IW. 

1036.  OKM 

STABILITY 

PERIOD:  lOSEC  ONE-  WAY  PROPAGATION  RUN 

MODE:  ORGIT 

TIME 

sa: 

ELL2TE 

SPECTRAL 

SCINTILLATION 

STEP 

(UT) 

AP 

LATITUDE 

;  LONGITUDE 

INDEX 

STRENGTH 

PHASE (RAD) 

AMPL(S4} 

1. 

2230 

40 

SO.  OON 

89.  lOW 

3.  9 

0.  0368 

1.  063 

0.  701 

2. 

2?31 

40 

79 

98.  OOW 

2.  9 

0.  0240 

1.  009 

0.  660 

3. 

2231 

40 

78.  87N 

109.  61W 

3.  9 

0.  0319 

0.  991 

0.  613 

4. 

2332 

40 

77.  71N 

119.  40U 

3.  9 

0.  0199 

0.  916 

0.  961 

S. 

3333 

40 

76.  2eN 

127.  3BW 

3.  9 

0.  0183 

0.  876 

0.  907 

6. 

2333 

40 

74.  64N 

133.  83U 

3.  9 

0.  0169 

0.  849 

0.  496 

7. 

2334 

40 

73.  87N 

139.  04W 

3.  9 

0.  0171 

0.  849 

0.  413 

8. 

2239 

40 

70.  f fN 

143.  30U 

3.  9 

0.  0193 

0.  900 

0.  389 

T. 

3239 

40 

6f .  04N 

146.  69U 

3.  9 

0.  0396 

1.  039 

0.  387 

10. 

2236 

40 

67  03N 

149.  83U 

3.  9 

0.  0662 

1.  671 

0.  467 

11. 

2337 

40 

64.  9BN 

192.  38U 

3.  9 

0.  3018 

3.  967 

0.  743 

12. 

2337 

40 

63.  fON 

194.  99U 

3.  9 

O.  3107 

2.  981 

0.  709 

13. 

2238 

40 

60.  79N 

196.  93U 

3.  9 

0.  0398 

1.  396 

0.  474 

14. 

223f 

40 

96.  66N 

198.  29U 

3.  9 

0.  0141 

0.  771 

0.  391 

19. 

233f 

40 

96.  91N 

199.  78U 

3.  9 

0.  0068 

0.  936 

0.  394 

14. 

3340 

40 

94.  36N 

161.  17W 

3.  9 

0.  0039 

0.  403 

0.  337 

17. 

3341 

40 

92.  19N 

162.  44W 

3.  9 

0.  0033 

0.  312 

0.  299 

IS. 

3341 

40 

90.  OIN 

163.  61U 

3.  9 

0.  0014 

0.  242 

0.  262 

If. 

3243 

40 

47.  6SN 

164.  69U 

2.  9 

0.  0008 

0.  181 

0.  217 

20. 

3243 

40 

49.  63N 

169.  69W 

2.9 

0.  0004 

0.  130 

0.  168 

Figure  8.  Sample  Output  (Orbit  Node) 


tennlnus  (approaching  geostationary  altitude,  for  Instance)  due  to  the  rudimentary  treatment  of 
the  tn-altu  drift  velocity,  of  the  Ionospheric  Irregularities  that  produce  radiowave 
scintillation.  This  deficiency,  which  was  most  severe  at  high  latitudes,  affected  the  calcu' 
lated  phase  scintillation  strength,  T,  through  the  dependence  of  T  on  the  effective  scan 
velocity,  V  ,  as  shown  In  Equation  (1)  of  Fremouw  and  Lansinger  (1981).  The  effective  scan 
velocity.  In  turn.  Is  a  function  of  the  scan  velocity.  V^,  of  the  propagation  line  of  sight  with 
respect  to  the  scintillation -producing  Irregularities  at  the  Ionospheric  penetration  point 
(Rino,  1979),  where  Is  a  function  of  the  satellite  orbital  velocity,  Vq,  and  the  (n-situ 
drift  velocity.  For  a  low-orbiting  satellite  (h  -  1000  km),  the  velocity  due  to  the  orbital 
motion  of  the  satellite  Is  normally  the  dominant  term  In  calculating  V^.  However,  for  slowly 
changing  (I.e.  high-orbit)  receiver-transmitter  geometries  (and,  as  we  will  show  later,  for 
some  rapidly  changing  geometries),  the  contribution  from  Is  comparable  to,  or  even  dominant 
over,  V^. 

The  drift  velocity  tnodel  used  by  HBMOO  was  a  very  rudimentary  one.  In  meters/sec.  It  was 
as  follows: 


Vdx-0 


Vdy  -50-15 


1  ♦  erf 


+  40(1  ♦  Xp) 


(15a) 

(15b) 


Vdj  "  0  (15c) 

where  erf  ■  error  function, 

Xj  -  Invariant  latitude, 

«  Invariant  latitude  of  high-latitude  scintillation  boundary, 

Kp  ■  3-hourly  planetary  magnetic  Index, 

and  X,  y,  and  z  denote  components  1n  the  geomagnetic  north,  east,  and  geocentric-nadir  direc¬ 
tions,  respectively.  Equations  (15)  describe  an  eastward  drift  of  50  m/sec  at  the  geomagnetic 
equator,  dropping  to  20  m/sec,  and  Increasing  with  geomagnetic  activity  at  latitudes  poleward 
of  with  a  maximum  velocity  of  20+80(1+Kp)  m/sec.  While  this  model  Is  fairly  representative 
of  at  low  geomagnetic  latitudes,  it  Is  largely  Inadequate  at  latitudes  poleward  of  the  high- 
latitude  plasma  trough. 

This  simple  model  produces  a  (magnetic)  eastward  drift  at  all  local  times  and  latitudes, 
whereas  the  existence  of  a  (nominally)  two-celled  plasma  convection  pattern  with  predominately 
anti-sunward  flow  across  the  polar  cap  with  predominately  sunward  flow  between  the  polar 
plasma-trough  wall  and  equatorward  of  roughly  75°  Invariant  has  been  accepted  for  at  least  the 
past  decade  (Cauffman  and  Gurnett,  1972;  Heells,  1982).  It  was  decided,  therefore,  to  survey 
the  various  two-cell  convection  models  that  have  appeared  In  the  literature  and  select  one  for 
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Implementation  In  Program  WBMOO. 

The  earliest  models  were  the  empirical  high-latitude  convection  potential  models  devel¬ 
oped  by  Heppner  (1972)  from  060-6  observations,  which  have  been  the  conceptual  ancestor  of  nrasi 
models  that  followed  (Wolf,  1974;  Volland,  1978;  Kawasaki,  1975;  Heppner,  1977;  and  Heells  et 

al,  1982).  Of  these  recent  models,  that  developed  by  Heells  et  at  (1982)  proved  to  be  the  most 
flexible,  both  phenomenologically  and  computationally,  and  ■wst  applicable  to  the  problem  'it 
hand. 

2-5.2  Convection  Potential  Model 

The  convection  potential  model  developed  by  Heells  et  al  (1982),  henceforth  denoted  the 
HLS  model.  Is  an  empirical  model  developed  for  use  In  an  Investigation  of  the  effects  of  the 
convection  flow  on  F-reglon  plasma  distributions.  It  Is  described  In  an  offset  Invariant  co¬ 
latitude,  (6),  rotated  magnetic  local  time  (0)  coordinate  system  In  which  all  local  time 
boundaries  are  assumed  to  be  along  lines  of  constant  6.  The  electrostatic  convection  poten¬ 
tial,  <|/,  Is  assumed  to  be  of  the  form 

4((^,6)  ■  G(0)F(4.,0)  (16) 

where  G(0)  and  F((j),0)  represent  the  latitudinal  and  magnetic  local  time  (W.T)  variations  of 
respectively. 

The  function,  G(3),  which  contains  the  main  latitudinal  variations  of  ij;,  Is  basically  a 
sin''©  form  with  r«2  In  the  polar  cap  and  r«-4  equatorward  of  the  potential  reversal  boundary. 
Two  transition  regions,  nne  on  either  side  of  the  reversal  boundary,  are  included  to  remove 
velocity  discontinuities  at  the  boundary.  The  exact  form  of  G(0)  Is  shown  In  Appendix  I.  In 
these  equations,  0^  is  the  co-latitude  of  the  reversal  boundary;  Oj,  and  8^  (0j  >  >  ©g)  are 
the  equatorward  and  poleward  boundaries  of  the  two  transition  regions;  and  A^,  Bj,  and  Bj 
are  scale  parameters,  which  are  determined  by  equating  the  functions  and  their  first  deriva¬ 
tives  at  6*0j  and  0*02,  respectively.  The  parameter  0^  is  an  additional  scale  parameter, 
denoted  as  a  polar  potential  phase  angle  by  HLS,  which  allows  for  a  non-zero  convection 
velocity  at  0«O°.  (The  determination  of  all  model  parameters  is  discussed  In  the  next 
section.)  Figure  9  shows  an  example  of  G(0)  for  0^*22*’,  0q»17®,  02*15’’,  0^*10°. 

Tne  function,  F((j),e),  which  models  the  MLT  variation  of  the  potential,  is  more  complex 
functionally  than  G(0),  but  basically  provides  for  constant  potential  along  ©"constant  except 
In  two  regions,  the  dayside  cusp  and  the  nightside  exit  region.  Appendix  I  gives  the  exact  form 
of  F((ti,0)  Implemented  in  WBMOO.  These  are  slightly  different  from  Equations  (5)  in  HLS  due  to  a 
slightly  different  Implementation  of  (MLT).  In  these  equations,  and  i{»g  denote  the  poten¬ 
tial  extrema  values  along  (J»»6”  and  <ii*18^,  respectively,  Is  the  MLT  location  of  the  center  of 
the  dayside  cusp,  (ti^  are  the  width  of  the  cusp  along  S«0q,  Is  the  MLT  location  of  the  night 
exit  region,  and  (>■  width  of  the  exit  region  along  Figure  10  shows  F(iti,0Q)  for 
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(0)9 


0  -35  kev,  *^—35  kev,  <>^-1200  MLT,  0*-  1^.  (>„»0200  MLT.  and  2^. 
m  e  fl  a  n  n 

Figure  11  is  an  example  of  the  polar-cap  potential  pattern  as  described  by  Equation  (16) 
using  the  model  parameters  employed  In  Figures  9  and  10.  Note  that  the  center  of  the  convection 
pattern  Is  offset  from  the  Invariant  pole.  This  Is  discussed  In  the  following  section  on  the 
selection  of  model  parameter  values. 

2-5.3  Model  Parameter  Definition 

One  of  the  most  attractive  features  of  the  HLS  model  Is  the  large  number  (16)  of  essen¬ 
tially  free  parameters  that  allow  the  user  to  configure  the  model  to  fit  observed  conditions 
within  the  two-cell  pattern  framework.  Table  1  Is  a  list  of  these  parameters  with  the  values, 
or  equations,  to  which  they  are  set  In  the  W6M0D  drift  model.  Also  Included  In  Table  1  are  the 
equations  used  to  calculate  three  environmentaT  parameters  (e,  Q^,  and  Kp)  from  Kp  If  they  are 
not  Input  to  the  model. 

The  philosophy  used  In  developing  the  values  or  functional  forms  shown  In  Table  1  was  to 
(1)  preserve  as  much  as  possible  of  the  essential  physics  of  the  convection  phenomenon  In  the 
choice  of  observables  and  model  paranibter  functions,  (2)  develop  expressions  that  would  drive 
the  model  based  on  routinely  available  observable  quantities,  and  (3)  allow  the 
entire  model  to  be  driven  by  a  single  geophysical  observable,  Kp,  If  no  other  data  are 
available.  This  last  criterion,  established  primarily  with  operationally  oriented  users  such 
as  AFGUC  In  mind,  led  to  development  of  the  expressions  shown  In  Table  1  for  e,  Qg,  and  Kp. 

Since  high-latitude  convection  Is  to  a  large  extent  "powered"  by  the  magnetosphere  and 
"braked"  by  the  Ionosphere,  It  seemed  reasonable  to  search  for  a  magnetospheric  parameter  to 
specify  the  "strength"  of  the  convection  (i.e.,  the  cross-polar  cap  potential  drop,  A\j;)  and  for 
auroral  Ionospheric  boundaries  to  specify  the  orientation  and  boundary  parameters  (I.e.,  H  q, 
4^1  9j.  “^d’  ***d’  ^n'  “^n^* 

The  most  widely  used  quantity  employed  to  specify  the  magnitude  of  magnetospheric  (and.  In 
turn,  high-latitude)  convection  Is  the  e  parameter  developed  by  Perrault  and  Akasofu  (1978)  to 
study  the  transfer  of  energy  from  the  solar  wind  to  the  earth's  magnetosphere.  Equation  (17)  1n 
Table  1  for  the  cross-nniar  ~C2p  pOvCMV  lal  drop  Is  taken  from  Relff  et  al  (1981),  who  studied 
several  parameters  to  employ  In  modeling  Ai|)  and  found  a  slightly  modified  version  of  e  to  be  the 
best.  In  this  equation,  e  Is  defined  as 

e  -  sln^  nT^-km/sec  (23) 

where  V  ■  streaming  velocity  of  the  solar  wind 

Sw  ^ 

^sw  *  solar  wind  8  angle 
B  ■  minimum  of 

®sw  *  solar  wind  magnetic  field  strength 
f  ■  amplication  factor  (■  7  for  use  In  Equation  (17)). 
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Figure  11.  Potential  pattern  displayed  In  Invariant  latitude-magnetic  local 
time  coordinates.  (Latitude  circles  are  every  10"  invariant; 
potential  contours  are  every  10  keV  with  a  maximum  of  +35  keV 
at  0600  HLT  and  a  minimum  of  -35  keV  at  1800  MLT.) 


TABLE  1 

POTENTIAL  MODEL  PARAMETERS 


Akasofu  solar  wind  paraaeter  c  -181. 7{K  +  3.08K  +  4.25)  (20)  Rieff  et  a1  (1981) 

^  Heppner  (1973) 
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The  selection  of  e  as  the  paraneter  on  Mhich  to  model  was  based  strongly  on  the  desire  to  use 
geophysically  meaningful  paraaMters  to  drive  the  model.  It  violates  to  some  extent,  however, 
the  desire  to  use  readily  (and  operationally)  available  parameters.  This  is  alleviated  in  the 
model  through  equation  (20)  in  Table  I,  in  which  c  is  calculated  as  a  function  of  Kp.  This 
relationship,  valid  only  in  this  application,  was  constructed  from  two  algorithms  for 

a/  -  0.93c  >  319  (Rieff  et  ol,  1981) 

All/  -  20  +  13  Kp  (Heppner,  1973). 

This  can  be  solved  for  c (Kp)  to  produce  equation  (20).  Note  that  Equation  (20)  is  strictly  a 
mathematical  artifact  which  allows  the  model  to  use  Equations  (17)  or  (24)  depending  on  the 
availability  of  e. 

Several  studies  have  been  made  of  the  partitioning  of  Ai|/  between  i|^  and  (Heppner,  1972; 
Heppner,  1973;  Rieff  et  at,  1981),  and  there  Is  evidence  of  systematic  asymmetries  between  |  | 
and  |i|<|^|,  which  can  possibly  be  tied  to  the  y-component  of  the  solar  wind.  By.  It  was  decided 
not  to  add  this  level  of  sophistication  to  the  model,  however,  since  no  quantitative  model  of 
the  variation  of  the  asymmetry  with  By  was  found.  Moreover,  B^  is  neither  readily  available  for 
operational  use  nor  can  it  be  simply  modeled  in  terms  of  some  available  parmeter  such  as  Kp,  as  c  is 
in  equation  (20).  After  reviewing  the  literature  on  this  point,  it  was  decided  that  4'g**'ji„*HAi|/ 
is  a  reasonable  default  partitioning  to  use  for  now. 

The  first  of  the  configuration  parameters  that  require  definition  is  the  location  (  H  q, 
<t>Q)  of  the  center  of  the  convection  pattern  in  invariant  latitude>MLT.  Although  Heelis  et  al 
(1980)  found  the  center  of  a  circle  fit  to  the  location  of  the  reversal  boundary  from  atmospher¬ 
ic  explorer  data  to  be  at  (4°,  22^40**),  it  was  decided  for  the  time  being  to  locate  the  center 
at  (5°  00^).  This  choice  was  made  for  the  following  reasons: 

1)  The  latitudinal  boundaries  (S2>^0>^l)  calculated  as  functions  of  various  auroral 
boundaries  based  on  auroral  oval  definitions  from  Holzworth  and  Meng  (1975)  and  Meng  et  al 
(1977),  who  locate  the  center  of  the  auroral  boundaries  at  roughly  (5°, 00**). 

2)  As  with  the  diifcision  to  Ignore  the  seasonal  variations  inherent  in  the  definitions  of 
MLT  (Section  2-1),  it  is  better  to  stay  simple  rather  than  to  add  unnecessary  (or  even 
fictitious)  variations  to  the  model. 

Several  different  ways  of  modeling  0^,  0q»  and  9j  were  explored  before  settling  on  the 
method  mentioned  earlier.  It  was  decided,  based  on  Figure  7  of  Heelis  at  al  (1980),  to  model 
the  poleward  edge  of  the  aurora  and  to  locate  the  reversal  boundary,  Oq,  2°  equatorward  of  that. 
The  expression  used  for  the  poleward  edge  of  the  aurora.  Equation  (18)  in  Table  1,  was  derived 
from  the  data  presented  Table  2  of  Holzworth  and  Meng(1975)  showing  the  poleward  edge  of  the 
statistical  auroral  oval  (Feldstein,  1963)  as  a  function  of  Q.  The  poleward  transition 
boundary,  02*  has  been  located  at  the  poleward  edge  of  the  oval,  as  both  Heelis  et  al  (1982)  and 
Volland  (1978)  placed  0^  roughly  2"  poleward  of  0q.  In  modeling  9^,  it  was  decided  to  use  this 
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parameter  to  control  the  drop  In  6(e)  below  as  a  function  of  the  equatorward  edge  of  the 
auroral  oval  (e^^).  Two  eiaplrlcal  equations  for  were  extracted  from  figures  In  Dandekar 
(1979) 

•  20.6  +  1.40  K  (25) 

ae  p 

-  22.6  +  1.04  Q  (26) 

from  which  Equation  (21)  In  Table  1  was  constructed  In  the  same  manner  (and  with  the  same 
caveats)  as  e(Kp)  was  derived.  Once  6^^  Is  calculated  from  either  Equation  (25)  or  (26),  6;^  Is 
determined  by  defining  It  as  the  value  needed  to  drop  6(6)  to  0.1  at  ^  If  there  were  no 
transition  zone.  This  leads  to  an  equation 

5ln  ej  ■  I  0.1(1  +  4  cot  Bj)**  I  ^  s1n(e^^  ♦  T)  (27) 

which  Is  solved  Iteratively  for  e^.  (It  should  be  noted  that  neither  transition  boundary  has 

any  physical  significance.  They  are  Included  only  to  avoid  discontinuous  behavior  of  Vq  near 

^e  reversal  boundary,  and  they  affect  the  model  potential  configuration  primarily  through 

slight  changes  In  the  equatorward  shape  of  6(e)). 

The  values  listed  In  Table  1  for  the  location  and  local-time  width  of  the  dayside  cusp  were 

taken  from  Heells  «t  al  (1976),  who  found  the  cusp  reversal  region  to  be  roughly  3>i  to  4  hours 

wide,  and  from  !4eng  (1981),  who  located  the  cusp  centered  at  roughly  11^30"VlLT.  No  evidence  was 

li<  iiiented  for  moving  the  cusp  center  MLT  ((f^)  with  magnetic  activity. 

The  location  of  the  night  exit  region,  however  does  apparently  move  toward  the  evening 

" jctor  with  Increasing  activity  (Z1  and  Nelson,  1982).  The  constants  In  Equation  (19)  In  Table 

were  developed  In  two  steps.  The  movement  of  i|>p  with  Kp,  the  mean  Kp  over  the  preceding  24 

4rs,  was  taken  from  Figure  4c  of  Z1  and  Nelson  (1982).  The  value  of  (jip  for  Kp  ■  Oo,  was 

determined  by  Iteratively  changing  until  the  model  produced  the  observed  location  of  as 

In  ''  Is  same  figure.  The  widths  of  the  exit  region,  are  currently  set  to  ±  2^  as  a 

conp.omise  to  the  two  examples  given  by  Heells  et  al  (1982)  In  their  Figure  6. 

For  completeness,  and  In  line  with  the  philosophy  stated  earlier,  a  relationship  was 

developed  between  Kp  and  Kp  (Equation  (22)).  This  was  developed  by  constructing  a  quadratic 

least-squares  fit  to  three  years  of  Kp  data  (1976-1978).  (Interestingly,  this  equation  Is 

nearly  Identical  to  one  developed  early  In  the  model  development  In  which  the  value  of  Kp  was 

assumed  to  tend  toward  3  as  one  went  back  24  hours  from  the  current  K.  value.  The  equation 

?  P 

from  this  simple  model  was  Kp  ■  1.1  +  0.60  Kp  -  0.016  Kp  .) 

The  power-law  parameters,  r^  and  r2,  were  set  to  values  suggested  by  Heells  et  al  (1982). 
The  value  for  e  suggested  (14°)  was  decreased  to  10°  to  avoid  undesirable  behavior  should  e  be 

C  w 

greater  than  Og. 
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Title  2  lists  the  vilues  for'Rp,  Q,  e,  6^,  0^,  •«<!  M  calculated  fron  the 

equations  presented  In  Table  1  for  Kp  values  0  through  9  .  Figures  12-14  are  examples  of  the 
convection  potential  pattern  produced  by  this  Implementation  of  the  HLS  model  for  Kp*2o,  3-f, 
and  60. 

The  functions  and  values  of  all  model  prameters  can,  and  most  likely  will,  be  changed  as 
the  model  Is  used  In  this  new  application.  As  will  be  discussed  below,  the  most  likely 
candidates  for  change  are  the  parameters  that  describe  the  location  and  width  of  the  night  exit 
region. 

2-5.4  Drift  Velocity  Calculation 

The  plasma  drift  velocity  Is  calculated  from 
-  E.xB 

Vd  •  (28) 

where  <l'(4,d)  Is  calculated  from  Equation  (16).  This  leads  to 


“^c  ^9®  r  39  r  sin 


(29) 


where  '^p  *  %  ■*■  hp  ■  6721  km.  From  Equation  (16) 


(30a) 


(30b) 


the  form  of  the  derivlatives  of  G(e)  and  F(^,0)  Implemented  In  Program  WBMOO  are  given  In 
Appendix  I.  The  coordinate  system  In  which  the  potential  Is  defined  has  Its  r-axis  aligned 
along  f,  so  Equation  (28)  becomes 


where 


Eq  -  -  7^  I  F(t,0)  +  G(9) 


E  ..-S4L_FJ[^ 

*^<1  rpsln  ■ 


(31) 


Figure  15  shows  calculated  from  Equation  (31)  plotted  In  the  (6,4)  potential -model  coor¬ 
dinate  system,  for  Kp«4o. 
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TABLE  2 
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Figure  13.  Same  as  Figure  12  for  Kp»3+. 
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The  final  step  in  the  calculation  is  to  transform  the  drift  velocity  calculated  in 
Equation  (31)  to  the  coordinate  system  used  in  WBMOD.  The  coordinate  system  of  Equation  (31) 

/N  ^  ^  A, 

has  '•-r  aligned  along  -B,  magnetic  south,  magnetic  east,  while  WBHOD  uses  magnetic 
north,  +y  magnetic  east,  and  +z  along  the  (geocentric)  nadir.  The  transformation  equations  are 


''dx  "  •''de  d  cose-  V^^sin  6  (33a) 

Vdz  =  ''de 

where  and  are  the  drift  velocity  components  in  the  potential  coordinate  system, 

Vdy,  ane  the  components  in  the  WBMOD  coordinate  system,  d  is  the  magnetic  dip  (inclination) 
angle  and  3  is  the  angle  between  x  and  3  given  by 


cos  6  = 


cos  Hq  -  cos  9p  COS  9j 
sin  ep  sin  6j 


(33) 


where  H  ^  =  5°  (from  Table  3),  0_  is  the  potential -model  0  coordinate,  and  0,  is  the  invariant 

0  P  •>  A 

co-latitude.  Figure  16  shows  the  V^|  plotted  in  Figure  15  after  the  transformation  shown  in 
Equations  (33).  This  figure  (and  others  that  follow)  is  ploted  such  that  the  0’  magnetic 
meridian  is  always  at  the  bottom  of  the  picture,  so  the  V^j  pattern  shown  is  rotated  according  to 
the  input  Greenwich  Mean  Time  (UT)  shown  on  the  figure. 

The  final  form  of  the  in-situ  drift  velocity  implemented  in  WBMOD  transitions  from  the 
drift  velocity  given  by  Equations  (33)  to  the  old  form  given  by  Equations  (15)  at  the  location 
of  the  sub-auroral  scintillation  boundary, (see  Section  3.2).  The  final  equations  are  then 


V 

V 

V 


dx 

dy 

dz 


+  p  V 
+  p  V 
+  p  V 


dxo 

dyo 

dzc 


(35a) 

(35b) 

(35c) 


where  the  primed  velocities  are  from  Equations  (33);  '^dzo  ®’"®  mid-latitude 

velocities  from  Equations  (15),  and  p  is  given  by 


where  Xj  is  invariant  latitude,  Xjj  is  the  scintillation  boundary  invariant  latitude,  and  Wj,  is 
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the  transition  width.  Figures  17  and  18  show  the  drift  velocities  calculated  from  Equations 
(35)  for  Kp  ■  30  at  1200UT  and  Kp  ■  60  at  1200UT  respectively. 

2-5.5  Effects  on  RMS  Phase  Calculations 

The  drift-velocity  model  described  above  recently  has  been  Implemented  In  WBMOD  (Rev  6B3), 
and  we  have  begun  to  assess  Its  effects  on  the  rms  phase  fluctuation,  ,  calculated  In  WBMOD. 
Results  to  date  are  tentative,  but  a  number  of  Interesting  manifestations  of  this  new  model 
have  turned  up. 

Figure  19,  a  plot  of  for  K  *3o  for  both  the  old  (6A3)  and  new  (6B3)  models.  Illustrates 
the  minimal  effect  expected  from  the  new  V^j  for  a  low-orbiting  satellite  case  (simulated  Poker 
Flat  night-overhead  Wideband  pass).  Figure  20,  however,  for  the  same  geometry  but  a  Kp  of  80, 
shows  an  unexpected  decrease  In  the  calculated  o^by  roughly  a  factor  of  two.  The  cause  of  this 
decrease  Is  the  addition  by  the  new  model  for  of  a  V^j^comparable  to  due  to  the  satellite 
motion,  which  drops,  by  roughly  45X,  the  effective  velocity,  Vg,  calculated  for  Kp"8o. 
(Recall,  from  Equation  15a,  that  V^^^-O  for  the  old  V^j  model).  The  source  of  this  and  the 
reason  that  little  change  occurred  for  Kp«3o  ,  can  be  seen  In  Figures  21  and  22.  These  figures 
show  the  track  of  the  penetration  point  plotted  over  the  drift  velocity  patterns  for  Kp“3o  and 
80  respectively.  In  the  Kp»3o  case,  the  track  lies  westward  of  the  exit  region;  the  majority  of 
the  drift  velocity  is  cross-track,  as  calculated  (roughly)  by  the  old  model.  For  Kp»8o  the 
track  lies  within  the  exit  region,  however,  and  the  drift  velocity  Is  almost  entirely  along- 
track,  effectively  reducing  the  Irregularity  scan  velocity. 

This  effect  Is  different  at  different  receiver  longitudes.  Figures  23  through  26  are  the 
same  as  19  through  23,  but  for  a  simulated  Goose  Bay  night-overhead  pass.  In  this  case,  the 
differences  In  f or  Kp»3o  are  more  noticeable  than  for  Poker  Flat,  and  the  decrease  In  at 
Kp=8o  Is  less.  As  can  be  seen  in  the  corresponding  pattern  plots,  the  cause  of  the  different 
behavior  Is  due  to  a  different  geometry  vis-a-vis  the  drift  velocity  pattern.  The  track  for  the 
Kp*3o  case  at  Goose  Bay  is  within  the  exit  region  and  has  an  appreciable  while  at  Kp»8o  the 
track  has  moved  closer  to  the  edge  of  the  exit  region,  so  that  the  along-track  drift  velocity  Is 

less  than  It  was  at  Poker  Flat  for  K„«8o. 

P 

In  another  test  of  the  model,  an  attempt  was  made  to  reproduce  the  results  of  a  study  (Basu 
et  al,  1982)  of  the  phase  scintillation  statistics  on  a  244-MHz  link  between  Goose  Bay  and 
FLEETSAT,  a  geostationary  satellite  at  (nominally)  lOO^W  longitude.  Figure  27  shows  a  plot  of 
median  diurnal  curves  of  observed  ajj,  for  Kp  >  3.5  and  average  model  for  Kp*6o.  The  dashed 
curves  are  for  Jan-Apr  1979  (upper  curve)  and  Aug-Nov  1979  (lower  curve)  from  Figures  3b  and  5b, 
respectively,  of  Basu  at  al  (1982),  and  the  solid  curve  Is  from  WBMOO  (Rev  6B3).  There  Is 
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Figure  19.  RMS  phase  for  a  simulated  Poker  Flat  night -overhead  Wideband  pass.  The 
solid  curve  Is  calculated  using  the  old  V.  model,  the  dashed  curve  using 
the  new  model.  ” 
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Figure  23.  Same  as  Figure  19  for  a  simulated  Goose  Bay  night-overhead  pass. 
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Figure  27.  RMS  phase  index  oa  for  a  Goose  Bay  -  FLEETSAT  geometry.  Solid  curve  is 
calculated  from  WBM0D-6B3  using  K  *60,  SSN»135.  Upper  dashed  curve  is 
from  Basu  eX  aJL  (1982)  Figure  3b,^loMer  dashed  curve  is  from  their 
Figure  5b. 


general  agreement  between  the  modeled  and  observed  but  some  disagreement  In  detail.  In 
particular,  the  peak  that  occurs  at  0900  UT  In  all  three  curves  is  caused  in  the  model  by  a  peak 
in  Vgff,  while  Basu  et  al  speculate  that  it  is  due  to  a  change  in  the  irregularity  strength,  C^, 
based  on  the  occurrence  of  a  similar  peak  in  S.  observations. 

In  summary,  preliminary  investigations  of  the  effects  of  the  new  model  have  produced 
several  unexpected  results.  If  correct,  they  should  be  observable  in  past,  such  as  Wideband, 
and  planned,  such  as  HILAT,  scintillation  experiments.  It  must  be  emphasized  that  the  most 
pronounced  of  these  effects  are  strongly  influenced  by  the  potential -model  parameters  that 
control  the  location  and  width  of  the  nightside  exit  region,  the  parameters  that  are  possibly 
the  most  inadequately  modeled  of  the  entire  set  in  the  present  model.  The  next  step  in  this 
process  should  be  a  careful  comparison  of  the  effects  predicted  by  the  new  model  with  the 
Wideband-Poker  Flat  data  base  to  adjust  the  model  parameters,  followed  by  comparirons  with 
other  data  sets,  such  as  the  Goose  Bay  -  FLEETSAT  observations,  for  verification. 
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SECTION  3 

WBMOO  REVISION  6B3 


3-1  OVERVIEW 

The  current  version  of  Program  WBMOD,  Revision  number  6B3  dated  31  January  1983,  Includes 
several  substantial  changes  to  the  Initial  version  described  by  Fremouw  and  Lansinger  (1981a}. 
The  major  changes  a:  e  as  follows: 

1.  Code  Structure  (6). 

a.  Geomagnetic  field.  The  basic  geomagnetic  field  model  (Subroutine  I6RF80)  has 
been  updated  with  the  International  Geomagnetic  Reference  Field  (IGRF)  1980  coefficients  as 
adopted  by  lAGA  In  1981.  Additional  routines  have  been  added  to  calculate  Invariant  latitude 
(Kluge,  1970)  and  local  magnetic  (dipole)  time  (Section  2-1). 

b.  Scan  velocity.  Two  modifications  were  made  In  the  calculation  of  the  llne-of- 
slght  scan  velocity  (Subroutine  VXYZ).  The  method  used  In  the  calculation  was  changed  tc  avoid 
an  artificial  contribution  due  to  rapidly  changing  magnetic  declination  at  high  latitudes,  and 
the  error  in  calculating  the  scan  velocity  discussed  In  Section  2-2  was  corrected. 

c.  Orbit  calculation.  The  circular-orbit  calculation  algorithms  discussed  in 
Section  2-3  were  Implemented. 

d.  User  Interface.  The  WBMOD  main  routine  and  subroutine  READIN  were  substantially 
modified  to  Improve  Interaction  with  the  user  and  to  make  changes  of  a  software-engineering 
nature  to  Improve  program  flow  and  code  transportability  to  other  computer  systems. 

2.  Model  form  (B).  The  only  change  in  model  form  from  the  original  WBMOD  is  the 
calculation  of  irregularity  drift  velocity  described  In  Section  2.5. 

3.  Model  constants  (3).  Three  model  constants  have  been  changed.  The  scintillation 
boundary  parameter,  (Section  3-3)  changed  from  71.0®  to  71.8®  invariant  latitude,  and  the 
height-integrated  strength  parameters  for  equatorial  (Cg)  and  middle  (C,„)  latitudes  were 
changed  from  2.3x10^  to  5.0x10^^  and  3.0x10^  to  1.9x10^^,  respectively. 

3-2  Structure  of  the  Code. 

figure  28  Is  a  flow  diagram  of  Program  WBMOD,  Revision  6B3.  Upon  Initiation  of  the 
program,  the  user  Is  asked  Interactively  for  Information  regarding  his  computational  scenario. 
The  requested  Information  Includes  parameters  of  the  user's  system,  such  as  operating  frequency 
and  the  longest  time  over  which  the  system's  mission  requires  phase  stability.  It  also 
includes  other  aspects  of  the  Intended  operation,  such  as  transmitter  and  receiver  location  and 
time  of  day,  plus  characterization  of  the  general  state  of  solar/terrestrial  disturbance  by 
means  of  sunspot  number  and  planetary  magnetic  activity  index,  Kp.  Finally,  the  user  specifies 
one  of  his  Input  quantities  as  the  Independent  varlab.e  (e.g.,  transmitter  location  or  time  of 


Figure  28.  Flow  chart  for  WBMOD  Rev  6B3. 
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day).  The  output  parameters  (e.g..  scintillation  Indices)  then  are  calculated  as  functions  of 
the  selected  Independent  variable. 

The  WBMOD  Driver  obtains  Information  froift  the  user  through  Subroutine  READ IN,  by  which  a 
mode  selection  also  Is  made  and  through  which  the  geometry  Is  initialized.  The  driver  then 
Increments  the  selected  Independent  varlable(s)  (Subroutine  STEP)  and  updates  the  geometry, 
calling  upon  the  Irregularity  model  and  a  propagation  theory  for  scintillation  calculations  at 
each  computational  point.  Among  the  parameters  that  the  user  may  choose  to  vary  are  the 
receiver  or  transmitter  latitude/ longitude  coordinates.  This  may  be  done  either  In  an  Incre¬ 
mental  but  static  manner  (mode  1)  or  In  an  orbital  mode  (mode  2)  In  which  the  scanning  motion  of 
the  line  of  sight  Is  taken  Into  account.  The  third  mode  consists  cf  varying  any  single 
Independent  variable.  (See  Section  3-4  for  a  more  detailed  discussion  of  user  Inputs.) 

The  output  parameters  -  value  of  the  Incremented  varlable(s),  p  (phase  spectral  Index),  T 
(phase  spectral  strength),  (RMS  phase),  and  S^  (Intensity  scintillation  Index)  -  are  written 
to  a  temporary  file  after  each  calculation  loop,  and  are  written  to  a  final  summary  file  (W80UT) 
by  Subroutine  PUTOUT  after  all  Increment  steps  are  completed  for  printing  or  plotting  by  means 
of  user-supplied  software.  When  the  chosen  parameter  range  is  satisfied,  the  user  may  start  a 
new  run  In  which  any  or  all  input  parameters  may  be  changed  without  a  total  reinitialization. 

Two  additional  output  files  can  be  generated  by  subroutines  BOUTM  and  BOUTS,  if  the  user 
desires,  by  configuring  subroutine  INIT  such  that  variables  LUBl  and  LUB2  are  Initialized  with 
FORTRAN  file  numbers  for  each  output  file.  Both  files  contain  the  penetration  point  latitude, 
longitude,  dip  latitude,  and  invariant  latitude  for  each  Increment  step.  In  addition,  the  file 
created  by  80UTM  contains  the  calculated  model  parameters  a  and  b  (along-  and  cross-field 
irregularity  axial  ratios),  (sheet  orientation  angle),  q  {in-situ  spectral  index),  C^L 
(height-integrated  strength),  and  hp  (height  of  the  phase  screen),  and  the  three  components  of 
In-situ  drift  velocity.  The  file  created  by  BOUTS  contains  auxllliary  and  scintillation 
parameters  p  (phase  spectral  Index),  T  (phase  spectral  strength),  G  (geometric  factor), 
(effective  scan  velocity),  o^,  and  S^. 

Table  3  contains  a  list  of  the  subroutines  and  functions  that  make  up  WBM00-6B3,  with  a 
brief  synopsis  of  the  function  of  each. 

3-3  MODEL  33 

Subroutine  MDLPRM  calculates,  or  calls  routines  to  calculate,  eight  model  parameters  used 
to  describe  ionospheric  irregularity  structure  for  the  calculation  of  the  scintillation  param¬ 
eters,  p,  T,  a^,  and  S^.  These  model  parameters  are  the  height,  l<p,  and  the  in-situ  drift 
velocity,  iTj,  of  the  irregularities;  the  outer  scale,  a;  the  height- integrated  spectral 
strength,  C^L;  the  In-aitu  spectral  Index,  q;  and  three  "shape"  parameters  describing  the 
three-dimensional  configuration  of  the  Irregularities,  a,  b,  and  6  .  The  philosophy  behind  the 
definition  and  form  of  these  model  parameters,  other  than  for  has  not  changed  from  the 
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Table  3 

Subroutines  and  Functions  In  Program  WBMOD 

Name  Synopsis 

WBMOO  Driver  Routine 

INIT  Initialize  :  .n  Routine  (user  specific) 

REAOIN  Calls  for  nputs  from  user  and  certain  subroutines 

AZNGCA  Finds  azimuth  and  grea^>c1rc1e  angle  between  points  1  and  2 
(Double  Precision) 

STCRB  Finds  clrr  lar  orbit  between  two  points  at  a  given  altitude 

HDLFRM  Establishes  a11  parameters  of  the  1onospher1c-1y  i^egularlty  model 
except  height-integrated  strength  and  in-aitu  drift  velocity 

FCSL  Computes  height-integrated  strength  of  Irregularities  from  empiri¬ 
cal  model 

ERF  Computes  error  function  (Double  Precision) 

IGRF80  Sets  up  arrays  for  calculation  of  the  International  Geomagnetic 

Reference  Field  (IGRF-1980}. 

SHELG  Calculates  geomagnetic  L  parameter 

STOER  Used  in  conjunction  with  SHELG 

FELOG  Calculates  magnetic-field  components 

GMLTC  Calculates  difference  between  local  mean  solar  time  and  local  mag¬ 
netic  (dipole)  time 

VDRIFT  Calculates  the  (n-aftu  drift  velocity  of  the  scintillation- 

producing  irregularities 

POTMOD  Calculates  nndel  parameters  fcr  the  high-latitude  convection 

potential  model  used  by  VORIF^ 

GFUNC  Calculates  the  latitudinal  variations  of  the  high-latitude  convec¬ 
tion  potential 

FFUNC  Calculates  the  local  magnetic  time  variations  of  the  high-lat.ltude 
convection  potential 

SOUTH  Builds  a  baseline  output  file  of  model  parameters  calculated  In 
MOLPRM 

ANGSD  Plane-geometry  routine 

COORD  Finds  point  2  given  point  1  and  azimuth  and  great-circle  angle 
between  them  (Double  Precision) 

VXYZ  Calculates  llne-of-sight  scan  velocity  (Double  Precision) 
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LATLOtt 

TRANS 

SCINT 

GAKMA 

HYP6E0 

CNUX 

6E0FAC 

OSRTN 

ROMINT 

FI 

FINDZ 

BOUTS 

STEP 

TSTEP 

PUTOUT 


Calculates  latitude  and  longitude  of  an  orbiting  satellite  (Double 
Precision) 

Calculates  coordinate  transformations  for  center-difference 
velocity  calculation  in  VXYZ  (Double  Precision) 

Calculates  scintillation  parmeters  using  phase-screen  scattering 
theory 

Ganma-function  routine 

Calculates  Gaussian  hypergeometric  function 

Calculates  normalization  factor  used  in  computation  of  intensity 
scintillation  index 

Calculates  static  and  dynamic  geometrical  factors  that  influence 
scintillation  strength 

Sets  up  integral  to  calculate  phase  variance,  for  finite  outer 
scale 

Modified  Romberg  quadrature  Integration  routine 

Computes  change  of  variable  for  efficient  integration  by  ROMINT 

Calculates  "reduced  height"  for  one-way  and  two-way  propagation 

Builds  a  baseline  output  file  of  scintillation  parameters 
calculated  in  SCINT 

Controls  the  incrementing  of  whatever  parameter  is  being  varied 
during  the  model  run 

Incraients  time  (includes  date  checks) 

Formats  and  controls  output 
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original  W6M0D.  (See  Section  III  of  Fremouw  and  Lansinger  1981a.)  The  new  calculation  Is 
discussed  In  Section  2-5.  Here  we  will  present  only  the  model  forms  for  hp,  a.  C^L,  q,  a,  b, 
and  6^  and  the  current  values  of  all  model  constants. 

1.  Irregularity  height. 


where 

hp  ■  Irregularity  height  at 

hg  ■  Irregularity  height  at  dip  equator 

^t  *  <^*e  ■  ^a)/2 

h,  ■  Irregularity  height  at  non-equator lal  latitudes 
X^  »  dip  latitude 

X^  ■  dip  latitude  of  b^-to-h^  transition 
■  half-width  of  transition. 

2.  Outer  scale. 

a  ■  1000.0  km  (effectively®). 

3.  Height-Integrated  strength. 
vtJT  -  E  +  M  +  H. 

a.  Equatorial  term,  (E). 


where 

Cg  >  equatorial  C^L  scale  parameter 
r 

er  >  Proportionality  constant  for  sunspot  number 
R  -  smoothed  sunspot  number 
>  ^  ■  dip  latitude 

X  g  >  dip  latitude  boundary  of  equatorial  "anomaly"  In  Ionospheric  total  electron 
content  and 


M 

Mhere 


H 

where 


Wg  ■  Half-width  of  equatorial-to-mid-latitude  transition 
Cgs  ■  proportionally  constant  for  season 
D  ■  day  of  year  (1-365) 

Aj  «  phase  delay  (days)  of  seasonal  variations 
*  geographic  latitude 

Cgg  *  proportionally  constant  for  geographic  latitude 
t  ■  local  geographic-meridian  time  (hours) 
tg  «  time  of  maximum  equatorial  scintillation 

■  post-maximum  temporal  half-width 
T_  -  T+  If  t  <  tg 

■  pre-maximum  temporal  half-width  If  t  >  tg. 
b.  Mid-latitude  term*  (H). 


+  Cgtcos 


exp 
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«  mid-latitude  C^L  scale  parameter 
■  variation  with  local  time 
t  *  local  time 
«  Invariant  latitude 

X^:  statistical  center  of  the  mid-latitude  region 
half-width  of  the  mid-latitude  region 

c.  High-latitude  term,  (H). 

-  Cut  cos  li  (t^  -  tj,^) 

“h  *  ^hb  b 

»  h1gh-1at1tude  C^L  scale  parameter 
>  proportionally  constant  for  sunspot  number 
R  *  smoothed  sunspot  number 
Xj  ■  Invariant  latitude 
X^  >  sub-auroral  scintillation  boundary 
W|^  *  width  of  transition  at  scintillation  boundary 
X^  *  nominal  Invariant  latitude  of  quiet-time  boundary 
Cjj  ■  rate  of  migration  of  ^  with  Kp  (deg  per  Kp  unit) 
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Kp  ■  three-hour  magnetic  Index  (0.0  to  9.0) 

Cbt  *  variation  of  with  local  magnetic  time 

*  local  magnetic  time  (LNT) 

W  *  phase  delay  (hours)  of  LMT  variation 
C^b  *  boundary-width  scale  parameter. 

4.  In-si tu  spectral  Index, 
q  -  1.5. 

5.  Field-aligned  axial  ratio. 

•  •  >0  ■  •.  [>  » 

where  a^  ■  equatorial  value  of  a 

*t  •  <®e  -  »a)/2 

a  *  non-equator la  I  value  of  a 

fl 

*  dip  latitude 

Xb  *  statistical  dip  latitude  of  a^-to-a^  transition 
Wb;  Half-width  of  transition. 

6.  Cross-field  axial  ratio. 

b  .  1  .  b,  [l  ♦  cos  (t,  -  t^)] 

where  bb  “  auroral  b  scale  factor. 

(See  high-latitade  C^L  term  for  other  variables.) 

7.  Orientation  angle  of  sheetllke  Irregularities  relative  to  L  shell. 

-  0.0. 

The  values  for  all  model  constants  In  Rev  6B3  of  Program  WBHOO  are  listed  In  Table  4.  Included 
In  this  table  are  the  constants  used  In  the  drift-velocity  calculation. 

3-4  USE  OF  THE  CODE 

Program  UBHOD  Is  structured  for  Interactive  application  from  a  user  terminal.  A  sample 
Interaction  for  Rev  6B3  Is  provided  In  Table  5,  In  which  system  queries  are  Indicated  In  lower 
case  and  user  responses  In  caps.  As  Indicated  In  the  table,  a  computation  session  begins  with  a 
request  by  the  code  for  a  label  by  which  the  run  output  Is  to  be  Identified.  The  label  may 
consist  of  any  alphanumeric  string  of  up  to  40  characters.  The  code  then  permits  the  user  to 
make  several  choices. 

First,  the  user  selects  either  one-way  (comounlcatlon  system)  or  two-way  (radar)  propaga¬ 
tion  and  then  the  reciprocal  of  the  low-frequency  cutoff  of  the  band  of  phase-fluctuation 
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TABLE  4 


HBHOO  Model  B3  Constants 

Irregularity  Height 

hj,  -  500  km  -  20.0®  dip 

Hq  »  350  km  Wf,  ■  3.0“  dip 


Height-Integrated  Strength 


Equatorial: 

5.0x10^^ 

Cer’ 

-  0.04 

10.0“  dip 

“e- 

10.0“  dip 

Ceg  -  15.0“ 

10.0  days 

22.5  hrs 

■  6.0  hrs  for  t 

T  «  T.  for  t  >  t 
-  +  e 

■  3.0  hrs  for  t 


Mid  latitude: 

«  1.9x10^^  -  0.33 

Wjjj  ■  15.0“  Invariant 
High  latitude: 

Cj,  «  4.3x10^^  *  0.0496 

Cbt  •  5.5  Cj,,  -  0.15 

tp,t  -  2.0  hrs 


»  32.5“  Invariant 

in 

Cj^  ■  1.5 

Xj  ■  71.8“  Invariant 


Axial  Katlos 

ag  «  30.0  a^  ■  8.0  •  0.75 

Drift  Velocity 

0.  ■  5.0“1nvariant  4  »  0.0  hrs  9,  »10.0“ 


WIMOO  RCV  M3  MQUCITI  A  CAICL  FOt  THIt  MN. 

POKER  FLAT  (N-OHI 

DO  YOU  WANT  ONK-UAV  1 1 )  DM  TWO-UAV  (31  AAOAADATtON? 

1 

INTOI  OUMATiON  (WCONM)  OVBI  UHICH  IvaTin  KOUlim  PHAW  OTAaiLITV 
(0.0  AON  OTDTCM  NOT  KNtITIWK  TO  PMAK  OCINTILLATIONI. 

10.0 

■NTtN  lONomnic  ouTn  ocaa^  (mo  or  tvpc  ‘NaosL*  rai  imicTtMLv  iminitv 

OOTAUCT  VALUC. 

NOOEL 

■Mm  tmwsuuMitTv  oRirr  mlocitv  (n/«)  rai  dkoharitk  loiTMi  iaot.  doum 
OR  typ«  ’nooii.'  niR  kr-ihrcnoknt  dkrault  wacuc. 

nOOEL 

RROVIDC  INITIM.  MN.UCO  ROR  ROLUMINO 

(IN  HMl,  DRClMAt.  NUNOIRi  OCCIMAL  HOURM.  DCCINAi,  CRRIil  OR  RM.  AO  ARRRBRRIATt). 
ORmATINO  nWaulNCV  (FRCOI: 

137.68 

RIAMKTARV  OUMAONtTIC  ACTIVITY  IM>I  (MR): 

3  0 

•nOOTIOO  ZURICH  RUNRROT  NUHRSR  (RRNI: 

75  0 

DAY  tr  TMK  YBAR  (DiWI: 

130.0 

HBRIDIAN  TII«  AT  RKCKIVOI  OR  RADAR  TAROKT  (TIM>: 

0.63 

lATITUM  OR  RSCtim  OR  TAROCT  (RtJkTt; 

63  1 


LONOITUDC  OR  RCCCIWfR  OR  TMOIT  (RljONI: 

-147.5 

ALTITVIOC  OR  RCCKIVCR  OR  TAROCT  (RALTt; 

0.2 

LATITUDC  OR  TRANSnirriR  OR  RADAR  COAT); 

80.3 

(.ONOITUDC  OR  TRANDHITTCR  OR  RADAR  (TLON): 

-03.7 

ALTITUOt  OR  TRANOHirm  OR  RADAR  (TM.T); 

1036  0 

RCLICT  CMAMINO  RARAHCTER 

(ONK  OR  TMK  NAMn  IN  RARRNR  AROVB  OR  TVRI  'RCRD*  OR  *TCRO* 
TO  STV  RtCCIVCR  OR  TRANRHITTER  ALONO  A  ORKAT-CIRCLK 
OR  ‘OROT*  TO  CAURC  HIOHCR  mniNAL  TO  RCRRORH  ORRITAL  OCAN 
AT  CONOTANT  ALTITVJOCI. 

QRBT 

RROVIN  rinal  value  ROR  the  FOLLOHINO: 

LATITUDE  OR  TRANDMITTER  OR  RADAR  (TLAT<: 

42.8 


LONOITUDE  or  TRANDHITTER  or  radar  (TLONI: 

-161.2 


MCIRY  MURDER  OF  CALCVLATION  ROINTO  (IN  (VUI: 

30 


dumhary  or  valuer 

RUN  CONOITIONR 

LAREL  .  ROKER  FLAT  (H-OHI 

t-3UAY  RH  OUR  OUTER  SCALC  VOZ  VOV  VOZ 

>  10.  O  MODEL  MODEL  DRIFT  VELOCITY 


INITIM. 

RREO: 

137.  M 
RLAT: 

AS.  100 


VALUES 

SON: 

73.  0 
RION: 
-147.  SOO 


DAY: 

IN. 

RALT: 

O.  300 


TIW: 

0.  AW 

TLAT: 

DO.  3M 


MR: 

3.0 

TLON: 

-03.  700 


TALT: 
lODA.  OOO 


CHMKOINO  RARAKETER  :  ORDT 

INITIM.  LAT  .'  SO  300  (JONO  ;  -03.  7M 

RINM.  LAT  :  43.  DM  LONO  :  -lAl.  3M 

4  OR  ROINTR  :  N. 

ARE  VALUER  CORRECT?  YER  >1  NO  •  3 

1 


DO  YOU  WISH  TO  start  AOAIN? 

I:  comrlete  rtartover 

3:  NO 

3:  ONLY  ALTERINO  ROME  VALUER 


2 


Table  5.  Semple  WBM0D-6B3  User  Interaction 


frequencies  to  which  his  system  is  susceptible.  Thereafter,  he  either  provides  values  for  the 
geomagnetic  east-west  outer  scale  of  the  in-situ  electron-density  spectrum  and  for  the  drift 
velocity  of  ionospheric  irregularities  or  elects  to  let  the  code  employ  default  values  for 
them.  For  high  latitudes,  the  default  values  for  drift  velocity  are  obtained  from  the  convec¬ 
tion  model  described  1o  Sections  2-5.2  through  2-5.4. 

Next,  the  code  requests  initial  values  for  the  potentially  variable  parameters  that 
describe  the  operating  scenario.  Once  the  initial  values  have  been  set  by  the  user,  the  code 
inquires  as  to  which  computational  mode  is  desired.  The  user  may  select  any  of  the  11 
parameters  Just  initialized  to  be  the  independent  variable  against  which  output  parameters  are 
to  be  tabulated.  Alternatively,  either  the  receiver  (radar  target)  or  transmitter  (radar)  may 
be  stepped  in  latitude  (north  positive)  and  longitude  (east  positive)  along  a  great  circle  by 
typing  in  RCRO  or  TCRD,  respectively,  instead  of  a  single  variable  name. 

In  eleven  of  the  foregoing  modes,  scintillation  parameters  are  calculated  for  each  incre¬ 
mented  value  of  the  Independent  variable(s)  without  introducing  a  line-of-sight  scan.  (That 
is,  scintillation  is  taken  to  arise  only  from  irregularity  drift.)  Finally,  the  ORBT  mode  may 
be  selected,  in  which  the  receiver  or  transmitter  (whichever  is  higher)  moves  along  a  constant- 
altitude  circular  orbit,  and  scintillation  results  from  a  combination  of  line-of-sight  scan  and 
irregularity  drift.  (In  most  low-orbiting  applications,  the  former  velocity  dominates.)  What¬ 
ever  mode  is  chosen,  the  code  now  asks  for  the  final  value(s)  of  the  changing  parameter(s)  and 
for  the  number  of  increments  desired  between  the  initial  and  final  values. 

A  sample  output,  (HBOUT),  corresponding  to  the  interaction  contained  in  Table  5,  is 
illustrated  in  Table  6.  Following  a  general  heading,  the  title  specified  by  the  user  is 
printed.  Thereafter,  his  input  parameters  are  identified,  followed  by  the  calculation  outputs. 
The  first  output  is  a  single  printing  of  the  power-law  spectral  index,  p,  of  phase  scintilla¬ 
tion.  Finally,  columns  containing  the  following  information  are  provided:  calculation  point 
number;  changing  paran»eter(s);  the  spectral  strength  parameter,  T,  for  phase  scintillation;  the 
phase  scintillation  index,  a^;  and  the  intensity  scintillation  index,  S^. 

Tables  7  and  8  illustrate  samples  of  the  model  and  scintillation  calculation  baseline 
output  files.  Both  outputs  include  the  penetration  point  geographic  latitude  (PLAT)  and 
longitude  (PLON),  dip  latitude  (DIPLAT),  and  invariant  latitude  (INVLAT).  The  baseline  model 
output  also  includes  the  calculated  values  for  a,  b,  6^,  q,  C^L,  Vp,  and  hp.  The  baseline 
auxilliary  and  scintillation  output  includes  the  calculated  values  for  p,  T,  6,  Vg,  Oij,,  and  S^. 


F-LAyn-FROCUCn  RAOIGMAW  KlNTtlXATION 


CAt.CUI.ATCD  ntOH  A  nOOCL  OCVCLORCB  CV  FNTSICAt.  DVMMtICC.  INC. 

ICLLCVUC.  UA  fMOR 


ucnoD  -  new  Ma 

THIC  Run  II  RORKR  FUAT  (M-OHI 
OtC-tlAV  PROFARATIQN 

RCOUIRCD  FHAK-ITARIUITV  DURATlCRt  -  to.  0  KC 
lONOVtCRIC  OUTCR  aCAUC:  CFRCCTIVELV  tlTINITC 
IRRCOULARITV  DRIFT  \«LOCITV  OCFAtR.T  HDOCU 


FRCI  «  137  «■  HHl 


RF  IMHI  •  a.  O  MN  •  75 


OAV  OF  VCAR  •  150. 


TINC  -  0.  AS  HOURI  LMT  AT  RCCCtWCR 


FOR  FIRST  DATA  POINT 

RCCCIUCR  COQRDINATCt  TRANIntTTXR  COORDUMTCS 

LAT  •  AS.  10  KC  uAT  •  SO  ao  K« 

LON  •  -147  SO  KC  LON  «  -05.  70  KC 

ALT  •  O.  aoo  AN  ALT  •  lOOA.  000  MM 


FOR  THIS  RUN.  TW  CHANOINO  PARAPKTCM  NCR| 

TRANRMITTCR  LATITUK  AND  LONCITUOC  SCM4NIN0  MJMC  KRIT 

POtCR-LAU  RPCCTRAL  INDCX  OF  PHAK  SCINTILLATION:  P  >  a.  50 


TLAT 

ILONO 

T 

8M1  PH4K(840) 

14 

1 

80.  200 

-89.  700 

0.  lAtAC-Ot 

0.  829 

0.  442 

2 

80.  033 

-91.  010 

0.  isaK-oi 

0  801 

0.  44!> 

3 

79.  744 

-94.  098 

0.  lAsat-oi 

0.  777 

0.  428 

4 

79.  434 

-100.  911 

0.  I34VC-01 

0.  794 

0.  409 

S 

79.  040 

-109.  412 

0.  tSTac-oi 

0.  733 

0.  990 

* 

78.  989 

-J09.  984 

0  ta(ME-OI 

0.  713 

0.  970 

7 

78.  077 

-113.  433 

0  tt4ac-oi 

0.  494 

0.  990 

a 

77.  933 

-114  941 

0.  iostc-di 

0.  477 

0.  929 

9 

74  929 

-120.  190 

0.  I03IC-0I 

0  443 

0  908 

10 

74.  399 

•123.  141 

0.  R«t7C-0a 

0.  449 

0-  487 

11 

79.  439 

-129.  837 

0  «ASAC-t>a 

0.  438 

0.  444 

12 

74.  993 

-138.  303 

0.  **««*-<« 

0  430 

0.  444 

12 

74.  243 

-130  941 

0.  lasAC-oa 

0.  434 

0.  424 

14 

73.  914 

-  J3.  433 

0  RiTac-oa 

0.  433 

0.  408 

19 

73.  747 

-134.  937 

0.  vaasc-oa 

0.  424 

0-  391 

14 

73.  009 

-134.  391 

0  laiae-oa 

0.429 

0.  374 

17 

71  339 

•137.  912 

0.  RSiK-oa 

0.  443 

0-  344 

18 

70.  443 

-139.  413 

0  lOTOC-Ol 

0.  473 

0.  399 

19 

49.  444 

-140.  804 

0.  taiAC-Oi 

0.  714 

0.  390 

30 

48  834 

-143.  103 

0  i4saE-oi 

0.  783 

0.  391 

31 

48  031 

-143  313 

0  ts7ac-ot 

0.  889 

U  340 

33 

47.  198 

-144.  449 

0  aAT|C-OI 

1.  041 

0.  381 

32 

44  348 

-149.  904 

0.  4ao4c-oi 

1.  331 

0.  423 

24 

49.  933 

-144.  904 

0.  7SAAK-OI 

1.  784 

0.  497 

39 

44.  493 

-147.  443 

0.  tAOtC  00 

a.  998 

0-  434 

34 

42.  844 

-148  330 

0.  asiAC  00 

4.  011 

0.  824 

37 

43.  997 

-149  149 

0.  SOSK  00 

4.  414 

0.  890 

28 

43  142 

-149.  944 

0  aaooc  oo 

3.  093 

0.  747 

39 

41.  389 

-190.  730 

0  SKTC-Ot 

1  927 

0  984 

20 

40.  434 

-191.  438 

0.  44taE-OI 

1.  344 

0.  482 

21 

99  940 

-193.  133 

0  assRGH>i 

1.  049 

0  439 

22 

98.  494 

-193.  777 

0  lATOC-OI 

0.  844 

0  389 

22 

97.  839 

-193.  403 

0.  laoac-oi 

0.  713 

0.  344 

24 

94.  992 

-154.  003 

0.  STsic-oa 

0  414 

0.  390 

29 

94.  080 

-194.  977 

0  AOMc-oa 

0.  939 

0.  338 

24 

99.  304 

-199.  139 

0  54aK-oa 

0.  478 

0.  339 

37 

94.  337 

-199.  440 

0  43ax-oa 

0.  428 

0.  331 

28 

93.  447 

-194.  173 

0  34541-08 

0  384 

0.  313 

39 

93  947 

-194.  644 

0  asa4c-oa 

0  349 

0-  309 

40 

91  489 

-197. 143 

0.  aas4C-oa 

0  310 

0.  294 

41 

90.  801 

-197  409 

0  tS47C-oa 

0  279 

0.  284 

42 

49.  914 

-198.  093 

0  i4aAe-oa 

0.  290 

0.  374 

43 

49.  030 

-198.  489 

0  1  tSK-oa 

0.  224 

0.  241 

44 

48.  143 

-158.  909 

0  54351-03 

0.  199 

0.  344 

45 

47  399 

-199.  314 

0.  74aAC-03 

0  177 

0.  230 

44 

44.  344 

-199.  711 

0.  375AC-03 

0  194 

0.  313 

4*» 

49.  479 

-140  097 

0.  44S7C-03 

0  138 

0.  199 

A8 

44.  989 

-140.  474 

0  a43AE-03 

0.  120 

0.  174 

49 

43.  493 

-140.  841 

0  aAI4C-03 

0.  109 

0.  197 

90 

43.  8r^ 

-141.  300 

0  I5S4C-03 

0.  091 

0.  139 

Table  6.  Sample  WBM00-6B3  Output 
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•••ASELINE  OUTPUT 


^8888g8g88888^Sg8888SS88S8888288888g8g8S8g8g88g8S38 

^!;8i!8888Sn888!i!iS:Ci88S8888i!!!t838888^:;88888S8;:!i:!:::: 

^§:-5SS52««22!is«??25S85fim335555i?5SS38j8 

‘88i8839S$380SHnS8i»;s88liSUi38S8a>H!8n???t? 

I  I  i  t  I  I  I  I  I  (  I  I  I  i  I  I  I  I  1,1  I  4  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  «  I  I  I  I  I 

n(«nc«nnMnrarin«inCiMOin(«n<^Nt«ei(«n(iii««(tfiiniifiiii«MnniPit«nfihiMni«it«rifi(v, 


p§§§§§§mi§mii§§3§i§§§ii§im§i§§i§§i§§§§i§§§§§ 

Sb6b66666Q666Qd6bdCt6ti666b66fi»6od66oo6b66666o6666o666 

3v»Anp'«nK«Ktr««*fii«»<««««i4iQ«*«g«no«»tto«*«»KnanNOir0tno99<^N««n 
i|4rT<w<rrlnn<1(lNtll«•«««oo>^V0^*^l'Q«»d««nN^l••••o^t>•KK<Kr»«vnr«rar•«• 
nncinnnnpifinnnnrihnridnfiioicinmfidoinwfifinnfinw-^MMMMM^****^*-*^ 

*!!!!!!!!!!!!!!^!!!!^!!!^!!!!!!!!!!!!!!!!!!!!!!!!!! 

^$s8a8§3&82ji2Si»SK22Si888SsSS2S22Ss8S322888S8ng3KS 

lKi«;!!!!;;;;!!SS;;;:SSt3338331l3SS333333333«3338};8S!igg88:iSS 

3d9Siii^l99i9****dninn<'bAbbb'«Nkkd00bbt>OQg-^«>«ninNnh«v\'‘AA«4 
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Table  7.  Sample  Model  Baseline  Output 


••■ASELSME  OUTPUT  (BCXNTI** 


Table  8.  Sample  Auxllllary  and  Scintillation  Baseline  Output 


SECTION  4 
RESEARCH  TOPICS 


4-1  INTRODUCTION 

The  scattering  theory  (Rino,  1979)  employed  In  WBMOD  is  formulated  for  a  simply  aniso¬ 
tropic  three-dimensional  power-law  spectral  characterization  of  the  scattering  Irregularities. 
By  "simply  anisotropic"  Is  meant  that  the  anisotropy  is  restricted  to  the  spectral  strength 
(and,  by  Inference  and  normalization,  to  the  Inner  and  outer  scales),  with  the  power-law 
spectral  index  taken  as  Isotropic.  The  theoretical  formulation  also  is  In  the  Infinite  outer- 
scale  limit. 

For  modeling  on  the  basis  of  such  a  characterization,  one  needs  six  spectral  parameters 
and  two  auxi Illary  quantities.  The  spectral  parameters  are  the  height-integrated  power  spec¬ 
tral  density,  C^L,  at  a  reference  wavenumber;  the  in-nitu  (one-dimensional)  spectral  index,  q; 
two  anisotropy-defining  axial  ratios,  a  and  b;  and  two  orientation  angles.  The  first  angle  is 
taken  to  be  the  local  magnetic  dip,  d,  on  the  assumption  that  the  Irregularities  are  elongated 
along  the  field  by  ambipolar  diffusion.  The  second,  6^,  Is  measured  from  the  magnetic  east  and 
is  invoked  to  describe  the  orientation  of  (sheetlike)  Irregularities  elongated  along  a  second 
axis. 

The  parameters,  a  and  b,  are  respectively  the  ratios  of  Irregularity  scale-size  along  the 
magnetic  field  and  along  a  second  elongation  axis  to  that  In  the  remaining  orthogonal  direc¬ 
tion.  The  two  auxllliary  parameters  needed  are  an  equivalent  phase-screen  height  (closely 
corresponding  to  the  centroid  height  of  the  Irregular  layer)  and  the  irregularity  drift  veloci¬ 
ty.  The  latter  is  needed  for  accurately  modeling  scintillation  on  geostationary  and  slowly 
moving  satellite  links.  It  was  addressed  in  Section  2-5. 

There  are  several  signal  parameters  observable  in  scintillation  measurements,  but  they  are 
by  no  means  related  one-to-one  to  the  foregoing  eight  Irregularity  parameters.  Nonetheless, 
relationships  do  exist,  and  a  relative  efficiency  can  be  achieved  by  applying  various  observ¬ 
ables  in  a  particular  order  to  e'.aluation  of  the  physical  parameters.  For  instance,  the  axial 
ratios,  a  and  b,  and  the  orientation  angle,  6^,  are  most  directly  deduced  from  interferometer 
measurements.  Approaches  that  do  not  use  interferometer  measurements  (Fremouw  and  Lansinger, 
1981a  and  1981b)  suffer  from  a  need  to  combine  multiple  observing  geometries  for  a  single 
result  and  from  uncertainty  over  variations  in  irregularity  strength. 

4-2  THE  PHASE-GRADIENT  RATIO 

General  analysis  of  interferometer  measurements  (Rino  and  Livingston,  1982)  is  quite  time- 
consuming.  For  application  to  modeling,  which  involves  large  volumes  of  data,  one  would  like  a 
more  time-efficient,  even  If  less  general,  technique.  In  our  research  related  to  WBMOD,  we 
explored  one  such  possibility,  the  theoretical  basis  for  which  Is  as  follows. 
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The  phase  structure  function,  D(^),  for  an  Interferometer  baseline  ^  is  defined  as 


D(6)  -<^[<t.(p+6)  -  tihyy  (36) 

and  is  related  to  the  normalized  phase  autocorrelation  function,  R(P),  as  follows  (Tatarski, 
1971); 


D(B)  «  2  [1  -  R(^)]  (37) 

2 

where  is  the  phase  variance.  For  a  power>1aw  phase  spectrum  with  spectral  index,  p,  and 
outer  scale,  a,  observed  in  a  processing  window  with  a  low-frequency  cutoff,  f^,  we  have  (Rino 
and  Fremouw,  1977;  Rino,  1979) 


R(8) 


J  uyi  '■>  *  fe V)  I"  * 


where  ■  a ^ 


(38) 


and  Kg  *  -y — 

and  where  is  the  effective  velocity  of  the  phase  pattern  across  the  observing  plane,  given  by 
Equation  (9)  in  Section  2-2. 

The  anisotropy  informaticn  we  seek  is  carried  in  the  argument  of  the  Bessel  function,  Jg, 
in  Equation  (38),  specifically  in  the  "effective  baseline"  given  by 


9,  ■  -  B6„a5,  *  /  (AC  -  4)‘‘  (3S) 

***'*’*®  ^'sx  *  ^  “  ®Z  ® 

Bsy  “  ^  tan  9  sin  <t>  (41) 


where  S  is  a  geometrical  factor  accounting  for  wavefront  sphericity;  A,  B,  and  C  are  given  in 
Rino  and  Fremouw  (1977);  eand  (|i  respectively  are  in  the  incidence  angle  and  magnetic  heading  of 
the  propagation  vector;  and  6  ,  S„,  and  6,,  are  the  components  of  the  interferometer's  physical 
baseline,  B.  Suppose  now  that  we  have  two  interferometers  oriented  perpendicular  and  parallel 
to  the  intersection  of  known  sheetlike  irregularities  with  a  plane  parallel  to  the  observing 
plane.  Calling  their  effective  baselines  B^and  B^^,  respectively,  we  can  define  the  "phase- 
gradient  ratio,"  r,  such  that 
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Our  desires  are,  first,  to  compute  and  B,,  for  different  Irregularity  axial  ratios,  a 
and  b,  and  observing  geometries,  which  dictate  A,  B,  and  C;  second,  to  evaluate  Equation  (38); 
and,  third,  to  identify  best  fits  of  the  ratio,  r,  to  pass-long  time  series  of  the  observed 
phase-gradient  ratio  (or  to  nultipass  averages  thereof).  Ue  would  like  the  p-dependenca  of  r 
to  be  sufficiently  weak  or  well-behaved  that  a  and  b  could  be  determined  from  the  behavior  of  r. 
This  would  give  us  a  means  for  evaluating  the  axial  ratios  that  Is  Independent  of  the  highly 
variable  Irregularity  strength,  but  much  faster  than  the  general  correlation  analysis  of 
Livingston  et  al  (1982).  We  must  first  Investigate  the  p-dependence  of  r  and  find  an  efficient 
means  for  evaluating  Equation  (38). 

Evaluating  the  Integrals,  whether  analytically  or  numerically.  Is  facilitated  by  making 
the  Infinite-outer-scale  approximation.  To  assess  the  effect  of  this  assumption  and  to  serve 
as  a  simple  Illustrative  exaiqple  of  the  behavior  of  Equation  (39),  we  first  evaluated  Equations 
(37)  through  (42)  for  the  following  Gaussian  autocorrelation  function: 

R(6)  =  exp(-6^/L2)  (43) 

where  Bg  Is  given  by  Equations  (39),  (40),  and  (41).  In  the  Infinite-outer-scale  approxima¬ 
tion,  we  have  obviously 

r  -  6i/B,|  .  (44) 

For  representative  model  parameters  and  actual  observing  geometries  and  baselines  employed  for 
Wideband  at  Poker  Flat,  we  found  the  general  Gaussian  case  to  yield  results  Identical  to  these 
obtained  from  Equation  (44). 

We  expect  the  Infinite-outer-scale  approximation  to  hold,  regardless  of  the 
autocorrelation-function  form,  so  long  as  Lq  Is  much  larger  than  the  physical  length  of  the 
interferometer  baseline,  and  we  have  employed  it  In  the  power-law  case.  Accordingly, 


so  we  have,  approximately. 


^  -ft) 

i(Bx.  p;  Kg) 

I(B||  »  P»  Kg) 

(46) 

where 

-  Jo(x)]  dx 

(47) 

xc 
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(48) 


and  . 


For  ccnputatlonal  efficiency,  mo  split  Integral  (47)  Into  two  ranges,  etrploying  respec¬ 
tively  small-  and  large-argument  approximations  to  the  Bessel  function.  Specifically,  we  used 

(49) 
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dx 
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(n+1)  -  x^(2n+l-p) 

4'’(nl)^(2n+l-p) 


(51) 
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“(p-*^)  [r2xCos  -  Tgy  Sin  I  (53) 

where  and  sre  the  real  and  Imaginary  parts  of  ^[-(p-Jj),  Ixjj],  and  and  are  the 
real  and  Imaginary  parts  of  r[-(p+*i),  Ix^^].  Ue  found  x^^  »  3.8  and  N  ■  5  to  give  quite 
satisfactory  accuracy  through  the  whole  range  of  Integration. 

We  coded  the  results  of  the  foregoing  analysis  In  the  model -development  version  of  WBMOD. 
Figures  29  and  30  contain  calculated  results  therefrom  to  Illustrate  pertinent  behaviors  of  the 
phase-gradient  ratio,  r.  The  left-hand  side  of  Figure  29  shows  that,  for  highly  anisotropic 
Irregularities  (bottom),  the  ratio  Is  sensitive  to  the  form  of  the  spatial  spectrum  (Gaussian 
vs.  power-law).  One  expects,  then,  that  the  ratio  also  may  be  sensitive  to  the  Index,  p,  of  a 
power-law  spectrum.  Indeed,  the  right-hand  side  of  Figure  29  shows  that  the  phase-gradient 
ratio  goes  essentially  as  the  effective-baseline  ratio  to  the  (p-l)/2  power. 

Figure  30  shows  graphically  the  sensitivity  of  the  phase-gradient  ratio  to  the  spectral 
Index  and  to  Irregularity  anisotropy.  It  contains  nine  plots  of  r  vs.  penetration-point 
position  during  the  pass  employed  In  Figure  29.  The  three  different  line  types  correspond  to 
three  different  values  of  the  cross-field  axial  ratio,  b,  for  fixed  field-aligned  axial  ratio, 
a,  as  Indicated  In  the  legend.  As  desired,  there  Is  quite  good  separation  of  the  curves  for 
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Figure  29.  Illustrating  (left)  the  sensitivity  of  the  phase-gradient  ratio,  r,  to  the  form 
of  the  spatial  spectrum  (Gaussian  vs.  power-law)  and  (right)  the  relative  impor¬ 
tance  of  the  two  factors  in  Equation  (46).  The  top  two  sets  of  curves  are  for 
Isotropic  irregularities  and  the  bottom  two  are  for  L-shell  aligned  (sheetHke) 
irregularities  with  axial  ratios  of  10:10:1. 


F  gure  30.  Illustrating  the  sensitivity  of  the  phase-gradient  ratio,  r,  to  the  cross-field 
axial  ratio,  b,  for  a  fixed  field-aligned  axial  ratio,  a,  of  10:1  and  to  the 
power-law  spectral  index,  p,  of  phase. 
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different  likely  values  of  b»  especially  for  the  mid-pass  geometry  In  which  the  line  of  sight  Is 
nearly  along  the  magnetic  field.  (The  peak  In  r  for  Isotropic  Irregularities  occurs  at  the 
geographic  zenith,  which  happens  to  be  close  to  the  geomagnetic  zenith  at  Poker  Flat.) 

Unfortunately,  dependence  of  r  on  p  Is  sufficiently  strong  as  to  render  b  ambiguous  as  a 
function  of  r.  Three  curves  are  shown  for  each  value  of  b,  employing  values  of  p  In  the  range 
observed  at  Poker  Flat.  For  r  >  1,  the  upper  curve  In  each  set  Is  for  p  »  2.9,  the  middle  curve 
Is  for  p  *  2.5,  and  the  bottom  one  Is  for  p  >  2.1.  (For  r  <  1,  the  order  of  the  curves 
reverses.)  Clearly,  unless  p  and  b  are  Independent  of  latitude  (which  they  probably  are  not), 
one  cannot  distinguish  a  situation,  say.  In  which  b  >  10  and  p  ■  2.1  from  one  In  which  b  >  5  and 
p  »  2.5. 

The  foregoing  shows  both  the  potential  utility  and  a  fundamental  limitation  to  employing 
the  phase-gradient  ratio  for  determination  of  axial  ratios.  It  appears  that  quantitative 
refinement  of  axial  ratios  from  a  large  data  population  may  be  feasible  using  the  phase- 
gradient  ratio,  but  that  underlying  variations  In  them  (e.g.,  latitudinal  dependence)  should  be 
modeled  with  guidance  from  more  general  Interferometer  analysis  (Livingston  et  ol,  1982). 
Moreover,  before  application  of  the  phase-gradient  ratio  for  this  purpose.  It  Is  necessary  to 
obtain  a  value  (or  a  functional  description)  of  the  phase  spectral  Index  Independently.  Phase 
spectral  behavior  Is  addressed  In  the  next  two  subsections. 

4-.'l  SPECTRAL- INDEX  BEHAVIOR  NEAR  THE  SUBAURORAL  ENHANCEMENT 

For  several  years,  we  have  known  that  the  values,  p,  of  the  phase  spectral  Index  observed 
In  Wideband  at  Poker  Flat  vary  with  geomagnetic  latitude  of  the  penetration  point,  at  least  at 
night.  Figure  31,  which  Is  adapted  from  Fremouw  and  Lansinger  (1979),  Illustrates  the  observed 
dependence.  The  Indicated  behavior  Is  potentially  of  Interest  for  modeling  because  of  the  one- 
to-one  relationship  between  p  and  the  powerlaw  Index,  q,  of  the  in-sltu  spatial  spectrum  of  the 
Irregularities  (q  ■  p  -  1  In  the  absence  of  diffractive  effects  on  the  phase  spectrum,  given  a 
thick  scattering  region).  We  did  not  pursue  It  Initially,  however,  because  of  a  concern  that 
the  observed  Increase  In  p  for  a  restricted  range  of  magnetic  latitude  might  be  an  artifact  of 
nonstationary  phase  statistics  due  to  rapidly  changing  observing  geometry  as  the  line  of  sight 
swept  through  the  region  of  geometrically  Imposed  phase-scintillation  enhancement  (Fremouw  et 
al,  1977;  Rino  et  al,  1978). 

Regardless  of  the  geophysical  reality,  or  lack  thereof,  of  the  latitudinal  dependence  of  p 
Illustrated  In  Figure  31,  that  dependence  Is  pertinent  to  Interpretation  of  phase-gradient 
measurements  In  terms  of  Irregularity  axial  ratios.  Moreover,  recent  Incoherent-scatter 
observations  of  large-scale  structures  In  the  F  layer,  apparently  restricted  to  the  latitude 
regime  near  Poker  Flat  (e.g.,  Vickrey,  1982,  Figure  1),  are  consistent  with  the  possibility 
that  Figure  31  reveals  a  true  geophysical  variation  In  q. 
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Figure  31.  Relationship  between  VHF  phase  spectral  index,  p,  and  geomagnetic  latitude 
(and  I  value)  at  the  1 ine-of-sight  penetration  point  in  the  F  layer. 
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In  view  of  the  foregoing,  we  conducted  a  preliminary  Investigation  of  the  behavior  of  p  In 
the  nighttime  Wideband  data  base  from  Poke*  Flat.  The  aggregate  behavior  and  an  Internally 
consistent  (hut  not  uniquely  proven)  Interpretation  of  It  are  described  In  this  subsection. 

Figure  31  shows  an  Increase  In  p,  at  night,  for  penetration-point  latitudes  corresponding 
to  magnetic  L  values  between  about  5  and  7.5.  Our  reason  for  doubting  the  geophysical 
significance  of  the  observation  stemmed  from  the  existence  of  sheetlike  Irregularities  1n  the 
nighttime  subauroral  Ionosphere.  The  line  of  sight  scans  through  grazing  Incidence  on  the  L 
shells,  along  which  such  sheets  ore  aligned,  at  L  ■  5.5.  Our  fear  was  that  the  scintillation 
enhancement  occurring  there  might  Invalidate  the  statlonary-statirtics-based  analysis  from 
which  the  measured  p  values  came. 

We  know  that,  as  a  dominant  configuration,  the  sheetlike  Irregularities  are  confined  to 
the  night  side  of  tiie  earth  (Fremouw  and  Lansinger,  1981b).  One  might  suppose,  therefore,  that 
a  test  of  the  geophysical  reality  of  the  observed  Increase  In  p  would  be  Its  presence  or  absence 
In  davtime  data;  absence  of  the  effect  In  the  daytime  would  Indicate  that  the  nighttime 
Increase  Is  a  geometrical  artifact  of  the  sheets.  This  reasoning  Is  oversimplified,  however, 
since  the  effect  may  well  be  associated  with  sheetlike  structures  but  still  have  geophysical 
significance. 

Nonetheless,  there  is  a  reliable  test  of  the  possible  effect  of  geometrically  Imposed 
signal-statistical  ncnstatlonarlty  on  our  observed  behavior  of  p:  to  look  for  a  singular 
feature  at  the  magnetic  zenith  In  our  daytime  data  base.  The  dominant  configuration  of  daytime 
Irregularities  In  the  vicinity  of  Poker  Flat  Is  that  of  axially  symmetric  "rods"  aligned  along 
the  magnetic  field  (Fremouw  and  Lansinger,  1981b).  For  such  Irregularities,  the  station's 
magnetic  zenith  represents  a  geometric  "singularity,®  which  Is  known  to  produce  an  enhancement 
In  phase-scintillation  strength  as  the  line  of  sight  sweeps  through  It.  If  the  same  rapid 
change  In  observing  geometry  does  not  produce  an  enhancement  in  p,  we  car  be  reasonably 
confident  that  our  nighttime  p  enhancement  Is  not  Imposed  by  a  geometrically  produced  signal- 
statistical  nonstatlonarlty. 

The  original  data  set  showing  a  broad  p  enhancement  (Figure  31)  contained  best-fit  spec¬ 
tral  Indices  to  15,953  phase  spectra  obtained  from  all  nighttime  VHF  Wideband  data  collected  at 
Poker  Flat  during  1976  through  1978,  except  for  exclusion  of  passes  during  which  the  local 
(College)  K  Index  was  zero.  Including  all  nighttime  data  (17,937  spectra  from  1976  through 
1979,  with  no  K-Index  threshhold)  produced  little  change  In  the  behavior  of  p.  The  latter 
result  Is  presented  In  Figure  32  (solid  curve)  In  terms  of  average  p  as  a  function  of  the 
Invariant  latitude  of  the  llne-of -sight  penetration  point  (at  350-km  altitude).  Shown  for 
comparison  (dotted  curve)  Is  a  plot  of  average  p  from  the  1,719  spectra  available  from  daytime 
passes  traversing  within  10°  of  the  magnetic  zenith,  which  occurs  at  a  penetration-point 
invariant  latitude  of  64.9°. 
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figure  32.  Average  VHP  phase  spectral  Index,  p,  as  a  function  of  Invariant 
latitude  of  the  line-of-sight  penetration  point  (at  350  kin)  for 
(a)  all  nighttime  Wideband  data  from  Poker  Flat  and  (b)  nearly 
overhead  daytime  passes. 


We  should  like  to  argue  either  (1)  that  lack  of  a  p  enhancement  In  the  dotted  curve  In 
Figure  32  attests  to  the  geophysical  significance  of  the  enhancement  clearly  observed  In  the 
solid  curve,  or  (2)  that  presence  of  an  enhancement  In  the  former  Indicates  that  the  nighttime 
enhancement  Is  a  geometrical  artifact.  The  results  shown  In  Figure  32,  however,  do  not  seem 
definitive.  There  are  more  statistical  fluctuations  In  the  dotted  curve,  owing  to  the  smaller 
number  of  data  points,  and  a  positive  fluctuation  does  occur  near  the  magnetic  zenith. 

In  an  attaapt  to  resolve  this  geometrical /geophysical  dilemma,  we  have  resorted  to  a  more 
directly  geometrical  ordering  of  the  data.  For  each  data  point,  we  have  calculated  the  angle 
between  the  line  of  sight  and  (1)  the  magnetic  L  shell  and  (2)  the  magnetic  meridian,  both 
calculated  at  a  height  of  3S0  km.  To  establish  whether  a  geometrical  "singularity"  due  to 
anisotropic  irregularities  produces  an  enhancement  In  p,  we  made  a  contour  plot  of  average  p 
values  on  a  grid  of  off-shell  and  off -meridian  angle.  The  question  at  hand  should  be  answered 
by  the  presence  or  absence  of  a  p  enhancement  at  the  origin  of  such  a  coordinate  system  for 
daytime  data. 

Figure  33  shows  such  a  contour  plot  for  all  the  6,647  daytime  spectra  available  within  the 
look-angle  window  employed.  There  Is  no  Identifiable,  statistically  outstanding  feature  at  the 
center  of  the  plot. 

The  grid  is  oriented  as  If  one  were  looking  down  on  a  map  of  Alaska.  That  Is,  positive 
off-meridian  angles  correspond  to  satellite  locations  to  the  (geomagnetic)  east  of  Poker  Flat, 
and  positive  off-shell  angles  correspond  to  northerly  satellite  positions.  Daytime  passes 
progressed  in  an  approximately  northwesterly  direction.  The  daytime  "overhead"  pass  corridor, 
from  which  data  for  the  dotted  curve  in  Figure  32  came,  starts  at  coordinates  of  about  +45°, 
-33°,  on  Figure  33,  progresses  slowly  thereon  through  the  lightly  dappled  region  (1.75  <ps  2.00) 
centered  near  +20°, -20°,  and  then  speeds  up  through  the  central  region  of  the  plot  (but  with  the 
corridor  center  passing  slightly  to  the  upper  right  of  0°,0°).  Tracing  such  a  trajectory,  one 
can  pick  out  the  features  that  rendered  the  dotted  curve  In  i^lgure  32  nondefinitive.  On  Figure 
33,  they  are  seen  to  have  no  special  geometrical  significance. 

Figure  34  contains  a  contour  plot  of  average  p  values  for  the  nighttime  Wideband  data. 
Unlike  Figure  33,  an  ordered  pattern  clearly  Is  present.  Surprisingly,  however,  it  is  ordered 
In  off-meridian  angle  as  well  as  In  off -shell  angle.  That  Is,  the  observed  enhancement  in  p 
does  not  appear  as  a  horizontal  baud,  as  one  might  have  expected  without  a  longitude-like 
separation  of  data  from  the  solid  curve  In  Figure  32.  Rather,  the  p  enhancement  occurs 
preferentially  away  from  the  station's  magnetic  meridian. 

The  Wideband  data  base  contains  about  twice  as  many  nighttime  data  points  as  It  does 
daytime  points.  In  order  to  assess  the  statistical  significance  of  differences  observed 
between  daytime  and  nighttime  behaviors  of  p,  the  nighttime  data  population  was  decimated  by 
using  only  every  other  chronologically  ordered  point  In  a  given  bin  prior  to  the  plotting  of 
Figure  34.  A  total  of  6,012  spectra  contributed  to  Figure  34  (as  compared  with  6,647  for  Figure 
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Figure  33.  Contour  plot  of  VHF  phase-spectral  index,  p,  for  all  daytime  Wideband  data  from 
Poker  Flit.  Gray  scale  appears  at  top.  Coordinate  system  refers  to  angle, 
calculat’d  at  350-km  altitude,  between  the  line  of  sight  and  the  local  (x) 
magnetic  meridian  and  (y)  magnetic  L  shell. 
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2),  and  the  numbers  of  points  per  bln  were  sii.<i1ar  through  the  Important  central  portions  of  the 
two  plots.  The  full  nighttime  daLa  population  {p  values  from  11,811  spectra)  was  employed  in 
Figure  35  In  order  to  refine  the  presentation  of  p  behavior. 

Figure  35  reveals,  even  more  clearly  than  Figure  34,  the  tendency  for  the  nighttime 
enhancement  In  p  to  avoid  the  magnetic  meridian.  Indeed,  this  Initially  puzzling  behavior  was 
observed  in  nighttime  corridor  data  sets  before  contour  plots  of  p  were  produced.  (Nighttime 
passes  progress  essentially  vertically  on  the  off -meridian/off -shell  coordinate  systen.) 

As  a  working  hypothesis,  the  following  interpretation  of  Figure  35  Is  offered.  Suppose 
that  the  "sheetllke"  Irregularities  (Rino  et  al,  19/8)  responsible  for  the  scintillation 
enhancement  (Fremouw  et  al,  19?7)  observed  by  means  of  Wideband  at  Poker  Flat  are  confined 
primarily  to  the  large-scale  end  of  the  scintillation-producing  spectrum.  That  is,  suppose 
that  the  smaller-scala  irregularities  are  more  nearly  axially  symmetric  than  are  the  large 
ones.  An  idealized  picture  would  be  ^small-scale  "rods"  Imbedded  In  large-scale  "sheets." 

The  scintillation  enhancement  ha«een  defined  mostly  In  terms  of  the  integral,  o|,  of 
the  phase  spectrum.  This  phase  variance  Is  thought  to  be  enhanced  by  quasi -coherent  addition 
of  phase  perturbations  as  a  radio  wave  propagates  along  an  elongated  axis  of  (sheetlike  nr 
rodlike)  Irregularities.  The  same  process  would  occur  for  the  rod/sheet  mix  postulated  in  the 
foregoing  paragraph.  For  such  a  mix,  however,  the  phase  spectrum  would  be  uniformly  enhanced 
only  near  the  magnetic  zenith,  while  there  would  be  preferantial  enhancenient  of  the  large-scale 
end  of  the  spectrum  at  large  off-meridian  angles  (end  small  cff-shell  angles).  For  the  latter 
geometry,  quasi -coherent  phase  addition  would  arise  In  the  (large-scale)  sheetllke  irregulari¬ 
ties  but  not  in  the  (small-scale)  radlike  Irregularities. 

The  foregoing  qualitative  reasoning  does  not  allow  us  totally  to  distinguish  geometrical 
from  geophysical  sources  of  the  p  enhancement.  We  assert  only  that  the  absence  of  a  p 
enhancement  at  0'’,0°  in  Figure  2  permits  us  to  accept  the  enhancement  shown  in  the  solid  curve 
of  Figure  32  as  "real"  —  i.e.,  as  other  than  an  artifact  due  to  geometrically  imposed  signal- 
statistlctl  nunstationarity.  We  cannot  yet  unequivocally  Interpret  the  observed  p  enhancement 
directly  as  a  latitudinally  ordered  q  enhancement.  Still,  there  Is  now  abundant  evidence  that 
the  nightside  F  layer  between  L  of  about  5  and  7  is  a  special  region  as  regards  plasma-density 
structure. 

A  likely  sequence  of  events  in  the  aforementioned  region  would  be  one  initiated  by  soft 
electron  precipitation  (Tanskanen  et  al,  1981)  setting  up  latitudinally  nonuniform  F-layer 
ionization  observed  by  Vickrey  (1982)  by  means  of  Incoherent  scatter  and  by  Rino  and  Owen 
(1980)  and  Leitinger  et  al  (1982)  as  TEC  enhancements.  This  L-shell -aligned  large-scale 
structure  would  form  a  "reservoir"  spectral  regime  from  which  smaller-scale  structure  would 
bifurcate  into  a  cascade  regime  by  means  of  convective  instability. 

As  a  new  element  in  the  picture,  the  p  behavior  reported  and  tentatively  interpreted  here 
may  indicate  a  trend  toward  cross-field  isotropy  In  the  cascade  regime.  Thus,  there  would  be  no 
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need  for  nonlinear  production  of  L-shell -aligned  structures  by  convective  processes,  which  do 
not  seem  to  be  a  dominant  result  of  NRL's  simulations  (Keskinen  and  Ossakow,  1982  and  refer¬ 
ences  therin).  One  may  envision  a  multistage  cascade  process  in  which  each  generation  produces 
structure  not  only  smaller  than  itself  but  also  oriented  perpendicular  to  its  own  dominant 
alignment  (parallel  to  its  own  dominant  gradient).  Many  generations  taken  together  would 
display  cross-field  isotropy,  but  the  largest  (reservoir-regime)  structures  would  remain 
anisotropic. 

For  propagation  calculations,  the  foregoing  picture  could  be  modeled  by  means  of  a  two- 
component  spectrum.  Signal-statistically,  the  situation  might  be  likened  to  the  two-component 
model  described  by  Fremouw  et  al  (1980),  in  which  the  large-scale  component  is  responsible  for 
large  phase  perturbations  and  geometrical -optics  focusing  and  defocusing,  while  the  small- 
scale  component  produces  diffractive  scatter.  Now,  however,  the  two  components  would  have  a 
phenomenological  basis  also,  distinguishable  not  only  in  scale-size  but  in  degree  of 
anisotropy. 

4  4  SPECTRAL  FEATURES 

The  spectral  index,  p,  behaving  as  described  in  Subsection  4-2,  is  obtained  (by  SRII)  from 
log-linear  fits  to  phase  spectra  between  0.5  and  10.0  Hz.  Designed  for  bulk  processing,  the 
fitting  procedure  presupposes  a  simple  (one-component)  power-law  spectrum.  Automated  fits  to 
more  complicated  spectra  could  produce  unexpected,  and  even  systematic,  behaviors  of  fit 
parameters  such  as  p.  Thus,  a  necessary  step  in  the  investigation  is  a  visual  inspection  of 
selected  spectra. 

Passes  for  initial  inspection  of  selected  spectra  have  been  chosen  in  the  following 
manner.  The  rectangular  boxes  on  Figure  35  represent  masks  used  for  convenient  selection  of 
pass  segments  for  comparison  of  the  disparate  behaviors  noted  near  and  away  from  the  local 
meridian.  In  terms  of  routine  20-second  processing  periods  for  which  we  have  detrend  data,  we 
found  17  points  from  six  passes  in  Region  1,  22  points  from  five  passes  in  Region  2,  and  28 

points  from  nine  passes  in  Region  3.  For  initial  inspection,  we  have  selected  periods  of  high  p 

from  Regions  1  and  3  and  periods  of  low  p  from  Region  2. 

Table  9  below  indicates,  in  terms  of  routine  20-second  periods,  that  the  points  obtained 

in  the  detrend-data  masking  procedure  are  consistent  with  the  behavior  of  p  in  the  overall 
(summary-data)  nighttime  population  from  Poker  Flat. 
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TABLE  9  SELECTED  DATA 


Region 

No.  of  Points 

Total 

p  >2.5 

p  <  2.0 

1 

17 

7 

3 

2 

22 

2 

13 

3 

2S 

8 

2 

We  selected  15  of  the  aforementioned  20  passes  for  detailed  spectral  analysis.  Including  a 
preponderance  showing  the  p-behavlor  Illustrated  In  Figure  35,  a  few  counter-examples,  and  one 
very  weak-sclntillatlon  pass  as  a  check  on  noise  effects.  For  each  of  the  15  passes,  we  plotted 
the  phase  time  series  for  five  or  six  minutes  Including  the  scintillation  enhancement.  The 
time-series  plots  were  used  to  hand-select  20-second  periods  showing  reasonable  statistical 
statlonarlty,  for  spectral  analysis.  Periods  totally  Inside  the  enhancement  region  and  totally 
outside  the  region  were  being  sought,  and  transition  periods  are  have  been  Identified. 

Figures  36  shows  some  of  the  spectra  computed  after  the  foregoing  selection.  The  cases 
were  chosen  as  exemplary  of  behavior  well  off  (a,  from  Region  3)  and  near  (b,  from  Region  2)  the 
magnetic  meridian.  The  former  (a.  Pass  61-33)  shows  an  ordered  increase  In  p  through  the  off- 
meridian  enhancement  region.  The  latter  (b.  Pass  52-46)  is  an  extreme  example  of  near-meridian 
behavior.  Not  only  does  p  remain  at  values  well  below  the  maximum  achieved  In  the  off -meridian 
case,  but  It  actually  decreases  substantially  In  the  strongest  portion  of  the  scintillation 
enhancement  (as  Identified  on  the  time-series  plot). 

The  spectra  Illustrated  In  Figures  36  were  obtained  using  a  Blackman-Harris  window.  The 
first  data  point  (0.1  Hz)  was  elevated  by  3  dB  to  offset  detrender  suppression,  and  a  five-bin 
centered  smoother  (modified  to  avoid  end  effects)  has  been  applied.  The  straight  lines 
superimposed  on  the  spectra  represent  log-linear  best  fits  In  the  spectral  band  routinely 
employed  for  that  purpose  In  Wideband  sunmary  processing.  The  spectra  shown  are  quite  well 
approximated  by  single-component  power  laws  In  that  band  and  above,  although  several  show  some 
enhancement  In  a  lower-frequency  (reservoir?)  band.  The  noise  floor  Is  at  about  -70  dB  on  the 
scale  Indicated.  The  p  values  derived  from  the  routine  fit  procedure  are  indicated  on  the 
plots. 

The  first  and  last  spectra  In  both  Figures  36  a  and  b  are  from  outside  the  geometrical 
enhancement,  and  they  exhibit  p  values  approximating  that  currently  employed  In  WBMOD  (2.5). 
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The  central  two  spectra  in  each  case  are  from  the  enhancements,  and  they  show  the  opposite 
spectral index  behaviors  described  above.  (Note  that  the  central  two  spectra  in  Figure  36b  are 
from  overlapped  time  periods  in  Pass  52-46,  which  displayed  a  narrower  enhancement  than  did 
Pass  61- 33.)  Qualitatively,  it  appears  that  the  off -meridian  enhancement  preferentially 
affected  large-scale  structures,  whereas  all  scales  were  enhanced  —  and  especially  small  ones 
in  the  present  example  —  near  the  magnetic  zenith.  (The  peak  enhancement  in  Pass  52-46 
occurred  as  the  line-of-sight  passed  within  about  two  degrees  of  field-alignment.) 

Reliable  interpretation  must  await  inspectiun  and  analysis  of  several  additional  cases, 
including  consideration  of  less  simple  spectra  than  those  shown  in  Figure  36.  Of  the  approxi¬ 
mately  one  hundreo  spectra  so  far  inspected,  about  three-quarters  display  some  general  spectral 
enhancement  above  the  routine  linear  fit  at  frequencies  below  0.5  Kz  (^  6  km  wavelength  in  the 
line-of-sight  scan  oirection).  The  second  most  comnon  departure  from  single-component  power- 
law  behavior  io  a  dowttward  break  at  the  high-frequency  end,  as  found  by  Rino  et  al  (1981)  in  in- 
si  tu  data  from  the  equator.  About  a  quarter  of  the  spectra  inspected  for  this  study  show  such  a 
break,  typically  at  about  ten  Hz  ( '  300  meters  wavelength  in  the  line-of-sight  scan  direction) 
or  slightly  lower.  One  of  the  most  obvious  examples  is  shown  in  Figure  37. 

In  the  case  shown,  which  is  from  Poker  Flat  Pass  12-19,  the  shallow-slope  band  happens  to 
coincide  with  the  band  used  for  routine  spectral  fitting.  In  other  cases,  such  coincidence 
does  not  occur,  and  the  best-fit  p  value  then  has  little  relevance.  A  single  power  law  is  still 
a  good  first  approximation  to  the  majority  of  spectra  inspected,  but  there  are  departures,  as 
noted . 

As  discussed  by  Rino  at  the  1982  DNA  Sunnier  Study,  it  is  tempting  to  interpret  spectra  such 
as  that  in  Figure  37  in  terms  of  "reservoir,  cascade,  and  diffusion"  regimes.  There  is  still  a 
good  deal  of  uncertainty  about  such  a  characterization,  however.  For  Instance,  a  3-dB  differ¬ 
ence  between  log-linear  best  fits  to  phase  spectra  in  bands  below  and  above  the  Fresnel 
frequency  could  be  explained  simply  in  terms  of  weak-scatter  diffraction  theory. 


FLUCTUATION  FREQUENCY  (Hz) 


Figure  37.  VHF  phase  spectrum  displaying  a  high-frequency  break  as  we11  as  a  low- 
frequency  band  that  Is  enhanced  with  respect  to  the  best  log-linear  fit 
between  0.5  and  10.0  Hz. 


SECTION  5 


CONCLUSION 

WBMOO  Rev  6B3  represents  the  state  of  the  art  of  descriptive  modeling  of  auroral -zone 
scintillation  for  systems-evaluation  use.  The  version  of  it  implemented  at  USAF  Global  Weather 
Central  contains  the  same  model,  the  code  being  configured  for  more  nearly  real-time  link 
analysis  on  the  basis  of  current  geophysical  input  conditions.  Neither  version  is  likely  to 
give  a  consistently  satisfactory  description  of  equatorial  scintillation. 

In  the  near  future,  WBMOO  Mill  be  iteratively  tested  against  equatorial  scintillation  data 
obtained  by  means  of  Wideband  at  Kwajaleln,  Marshall  Islands;  and  Ancon,  Peru.  This  work  will 
be  guided  by  the  comparative  morphology  of  Livingston  (1978)  and  the  physical  insights  of 
Oachev  and  Walker  (1983)  and/or  Tsunoda  (1983). 

The  drift- velocity  model  described  in  this  report  requires  additional  calibration, 
particularly  of  its  behavior  in  the  night  exit  region.  The  existing  Wideband-Poker  Flat  data 
base  will  be  used  in  this  calibration,  as  well  as  data  from  the  various  auroral  incoherent 
scatter  radars  (Foster,  1983),  and  from  the  HILAT  experiment  when  its  data  become  available. 

Meanwhile,  research  will  continue  into  the  three-dimensional  spectral  characteristics  of 
scintillation-producing  Irregularities  in  the  auroral  F  layer.  Current  information  on  the 
latitudinal  dependence  of  irregularity  axial  ratios  (Livingston  et  at,  1982)  will  be  incor¬ 
porated  into  WBMOO.  Beyond  that,  an  attempt  will  be  made  to  ascertain  whether  cross-field 
anisotropy  may  be  restricted  in  scale  size  as  well  as  in  geomagnetic  latitude  and  time. 

The  three-dimensional  spatial  spectrum  of  plasma  structures  is  expected  to  hold  important 
clues  about  their  generation,  evolution,  and  decay.  Thus,  it  will  continue  to  be  explored  with 
Wideband  data.  Moreover,  the  ability  of  scintillation  measurements  to  yield  such  informaton  is 
expected  to  constitute  a  major  contribution  to  the  beacon  experiment  in  the  forthcoming  HILAT 
Satellite  Program. 
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